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ABSTRACT 

We investigate the evolution of the H/3+[OIII] and [Oil] luminosity functions from z ^ 
0.8 to ^ 5 in four redshift slices per emission line using data from the High-z Emission 
Line Survey (HiZELS). This is the first time that the H/3+[OIII] and [Oil] luminosity 
functions have been studied at these redshifts in a self-consistent analysis. This is 
also the largest sample of [Oil] and H/3-|-[OIII] emitters (3475 and 3298 emitters, 
respectively) in this redshift range, with large co-moving volumes ~ 1 x 10® Mpc“® in 
two independent volumes (COSMOS and UDS), greatly reducing the effects of cosmic 
variance. The emitters were selected by a combination of photometric redshift and 
color-color selections, as well as spectroscopic follow-up, including recent spectroscopic 
observations using DEIMOS and MOSFIRE on the Keck Telescopes and FMOS on 
Subaru. We find a strong increase in L* and a decrease in for both H/3-|-[OIII] and 
[Oil] emitters. We derive the [Oil] star-format ion history of the Universe since z ~ 5 
and find that the cosmic SERB rises from z ~ 5 to ~ 3 and then drops towards 
z ~ 0. We also find that our star-formation history is able to reproduce the evolution 
of the stellar mass density up to z ~ 5 based only on a single tracer of star-format ion. 
When comparing the H/3-|-[OIII] SFRDs to the [Oil] and Ha SERB measurements 
in the literature, we find that there is a remarkable agreement, suggesting that the 
H/?-|-[OHI] sample is dominated by star-forming galaxies at high-z rather than AGNs. 

Key words: galaxies: evolution, galaxies: high-redshift, galaxies: luminosity function, 
mass function, cosmology: observations 


1 INTRODUCTION 


Our understanding of the mass assembly and star-formation 
processes of the Universe has improved greatly over the past 


few decades (for in-depth reviews, see Kennicutt & Evans 
2012 and Madau fc Dickinson|2014 1. We currently have ev¬ 


idence to show that the cosmic star-formation rate (SFR) 
peaked at z > 1 and that about half of the current stellar 
mass density had been assembled by that time (e.g, Lilly 


|et al.|p'996| [Hopkins fc Beacom||200^ [Karim et al.|[^ll[ 
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Sobral et al.||2013 l. However, many open questions remain. 


How fast did the star-formation rate density (SERB) drop 
at z > 2? How has the population of star-forming galaxies 
changed over cosmic time? How does the evolution depend 
on the enviroment over cosmic time? To answer these ques¬ 
tions, it is imperative that we use samples of star-forming 
galaxies that are low in contaminants and are well-defined 
in terms of selection methodology. 


There are many different star-formation indicators and 
calibrations in the literature. Each indicator traces the star- 
formation activity in galaxies independently and with differ¬ 
ent timescales. The ultraviolet (UV) light from bright, young 
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stars with masses > 5 Mq traces the bulk of the young pop¬ 
ulation with time scales of ~ 100 Myr. UV light of stars 
with masses > 10 Mq ionize the gas along the line-of-sight, 
resulting in the absorption and then re-emission of photons 
seen as nebular (e.g., Lyman, Balmer, and Paschen series 
of the Hydrogen Atom) and forbidden emission lines (e.g., 
[OIII] & [Oil]). The lifetimes for stars capable of ionizing the 
surrounding gas and dust to form the nebular and forbid¬ 
den emission lines are on the scale of ~ 10 Myr, allowing for 
the measurement of the instantaneous star-formation rate. 
Other indicators include far-infrared emission coming from 
the heating of dust shrouding the hot, UV bright, young 
stars and the synchrotron emission in the radio coming from 
accelerated electrons in supernovae. For an in-depth review 
of the various indicators and calibrations, we refer the reader 


to reviews in the literature (e.g., Kennicutt 1998 Calzetti 


Despite the different indicators that exist, one can not 
say that only one of the indicators is the “holy grail” of 
measuring the SFR of star-forming galaxies. However, us¬ 
ing different tracers to map out the evolution of the cosmic 
SFR history is not the best solution either. This is because 
evolutionary studies based on samples selected in different 
ways and using different indicators/calibrations at different 
redshifts will be susceptible to complicated, strong biases 
and selection effects, which results in a significant scatter 
when combining all of them to probe the evolution of the 
cosmic SFR. Another issue is that most studies don’t probe 
sufficiently large volumes to overcome the effects of cosmic 
variance. Furthermore, the effects of correcting for dust ex¬ 
tinction, especially for UV and optical studies, can result 
in large systematic uncertainties. One requires an indicator 
that can be used to probe from the \ow-z to the high-z Uni¬ 
verse using a robust and consistent methodology to reduce 
the effects and biases that come from making assumptions 
and differing selection techniques. 

Emission lines observed using narrow-band imaging 
techniques can provide an accurate and reliable sample of 


star-forming galaxies (e.g., 

Bunker et al. 

1995 

Fujita et al. 

2003 

Glazebrook et al. 

|2004| Ly et al. 

2007) 

Villar et al. 

2008 

Geach et al. 2008 

Sobral et al. 20131. The method- 


ology utilizes two different images of the same field: one 
being from a broad-band filter and the other being from a 
corresponding narrow-band Hlter. The narrow-band image 
is dominated by emission-line galaxies and the continuum, 
while the broad-band image is dominated by the continuum 
with a small contribution from the emission-line. When the 
two are subtracted, the result is the removal of the contin¬ 
uum and an image of galaxies with emission-lines. The ad¬ 
vantage of narrow-band imaging surveys is that they allow 
for the selection of emitters with a clean selection function 
by emission line flux and within a known narrow redshift 
range. This is because the Hlter width is quite narrow, such 
that any source brighter than expected from its broad-band 
magnitude is an emitter. Emission-line surveys via grism 
spectroscopy on the HST (e.g., Colbert et aI.|20T3 l are also 
great accompaniments to narrow-band studies, as they are 
area-limited (area of the grism) while emission-line surveys 
are redshift-limited. 

Most narrow-band surveys have focused on Ha (e.g.. 


formation indicator which is well-calibrated in the local uni¬ 
verse and is only mildly affected by dust attenuation. The 
latest results of the High-Emission Line Survey (HiZELS, 
Geach et al.||20d8 Sobral et al.||20()9 2012 20131 have ro¬ 
bustly traced the evolution of the cosmic SFR up to a ~ 2. 
This is the maximum redshift that Ha surveys can probe 
from the ground, as at higher redshifts Ha falls into the 
mid-IR and is blocked by water vapor and carbon diox¬ 
ide in the atmosphere. To probe to higher-a using the same 
narrow-band technique would require another emission line. 
The other major emission lines associated with star-forming 
galaxies are H/34861, [OIH]4959, [OIHJ5007, and [OII]3727 
which can be probed up to a ~ 3 for H/3-|-[OHI]and up to 
z ~ 5 for [Oil]. 

In the past decade, several H/1, [OIII], and [OH] stud¬ 
ies have been carried out (e.g., 


Hammer et al. 1997 


Hogg 


et al. 

1998 

Gallego et al. [2002 Hicks et al.| 

2002 

Teplitz 

et al. 

2003 

Ly et al. |20071 Takahashi et al.| 

2007 

Bayliss 

et al. 

2011 

2012 

Sobral et al. 2012| Giardullo et a 

. 

Drake et al 

2013 

1, the majority of which had small sam- 


Tresse et al 

120021 IFujita et al. 12003 

|Pascual|2005[|Ly et al.| 

to 

0 

0 

2011 

Sobral et al.||2009 |MT3 

1 as it is a reliable star- 


pie sizes and, hence, suffered from cosmic variance biases. 
The majority observed up to a ~ 1, while only the works of 
[Bayliss et al. (20111 and Bayliss et al. (20121 measured the 
[OH] SFR densities at z ~ 2 and 2 : ~ 4.6, respectively. Both 
these works used small samples, with the 2 : ~ 4.6 measure¬ 
ment having a sample size of only 3 [Oil] emitters, greatly 
limiting any conclusion. 

This paper presents, for the first time, the luminosity 
functions of H/?-|-[OHlQand [OH] emitters up to 2 : ~ 5 using 
the reliable selection techniques of |Sobral et al. (2009 2012 
|2013P on the combined COSMOS and UDS narrow-band 
publicly available catalog^^f HiZELS. The sample probes 
comoving volumes of up to ~ 1 x 10® Mpc®, which greatly 
reduces the effects of cosmic variance. This is also the largest 
sample of H/?-|-[OHI] and [OH] emission line galaxies up to 
2 ; ~ 5 to date in the literature and is used to effectively and 
robustly probe the evolution of the cosmic SFR density. 

This paper is organized as follows: in Section 2, we 
outline the photometric and spectroscopic data sets that 
we use and the methodology utilized to effectively select 
H/3-|-[OIII] and [OH] emitters. In Section 3, we outline our 
volume calculations and the completeness and filter profile 
corrections. Also in Section 3 are the results of the luminos¬ 
ity functions. In Section 4, we present the results and discuss 
our cosmic star-formation rate densities and the evolution 
of the stellar mass density based on [Oil] emitters. Section 
5 outlines the conclusion of this work and is followed by ap¬ 
pendix A, which presents our color-color selection criteria, 
and appendix B, which presents our binned luminosity func¬ 
tion data points, and appendix C, which presents our SFR 
density compilation. 

Throughout this paper, we assume ACDM cosmology 


^ Because the H/3 and [OIII] emission lines are close to each other, 
photo -2 and color-color selections can not distinguish between 
them. The best way to fully differentiate the two is via spec¬ 
troscopy. Based on line ratio studies, we can argue that most of 
the emitters will be [OIII] emitters, but to ensure we are not bi¬ 
asing our measurements based on such assumptions, we present 
the results as the combined measurement of H/?-|-[OIII]. 

^ The narrow-band catalogs are available on VizieR and are from 


Sobral et al. (2013l. 
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with Ho = 70 km Mpc“^, = 0.3, and f^A = 0.7 
with all magnitudes presented as AB and the initial mass 
function is assumed to be a Salpeter IMF. 


2 DATA 

2.1 Selection Catalogs 


Our data consist of narrow-band and broad-band photomet¬ 
ric data of the UDS ( [Lawrence et al.||2007| ) and COSMOS 
I Scoville et al.|2007 1 fields and spectroscopic follow-ups that 
are described in sectio n |2.4| The cat alogs that are described 
below are taken from jSobral et al.j ( |2013| ) and are publicly 
available. 

The narrow-band catalogs are from the High -2 Emis¬ 
sion Line Survey (HiZELS, [Geach et al.|[2008||Sobral et al.j 
2009 2012 20131. This project utilizes the narrow-band J, 


H, and K filters of the Wide Field CAMera (WFCAM) on 
the United Kingdom Infrared Telescope (UKIRT) and the 
NB921 filter of the Suprime-Cam on the Subaru Telescope. 
Previous uses of this data focused primarily on Ha emit- 


ting galaxies and their properties (e.g.. 

Cam et al. 

2010 

Geach et al.|2012| [Sobral et al.|2010l|2011| 2013 

Stott et al.| 

2013a|b 

Swinbank et al.[[2012a|b 

Darvish et al. 

2014 

1, but 


HiZELS is able to pickup more than just this emission line as 
it can detect any line above the flux limit. In this paper, we 
focus on the H/3-|-[OIH] and [OH] emitting galaxies found by 
HiZELS. The broad-band catalogs are from the Cosmolog¬ 
ical Evolution Survey (COSMOS; [Scoville et aL||2007 l and 


the DR8 release of the Subaru-AMM-UKIDSS Ultra Deep 
Survey (UDS; Lawrence et al.|2007[) catalog]^ 


2.2 Narrow-Band Selection of Potential Emitters 


In this subsection, we will review the methodology of select¬ 


ing potential emitters followed by Sobral et al. (20131. We 


refer the reader to this paper for a detailed overview. 

Potential emitters were selected by their color excess in 
terms of the parameter E (Bunker et al.[[T995 Sobral et al. 


2012), which quantifies the significance of the excess of a 


source with respect to the random scatter expected for a 
source to have a color excess > 0 in terms of their narrow- 
band magnitudes: 

1 _ ]^q-0.4{BB-NB) 


E = 


10-0'4(zp-nb) V^rip(4B - <^iB) 


( 1 ) 


where BB and NB are the broad-band and narrow-band 
magnitudes, respectively, rap is the aperture radius in pixels, 
and (tbb and unb are the rms per pixel for the broad-band 


and narrow-band, respectively (see e.g. Sobral et al. 2012 

[M^ . 

Emitters are selected on the basis that they have E > 3 
and a rest-frame equivalent width EWo > 25 A. The second 
condition ensures that we select sources with a significant 


^ We refer the reader to the respective papers for further de¬ 
tails on the creation of these catalogs. We also refer the reader 
to the UKIDSS and COSMOS websites for further information of 
the multi-wavelength photometric and spectroscopic data sets. 
(COSMOS: http://irsa.ipac.caltech.edu/data/COSMOS, UDS: 
http://www.ukidss.org) 
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Figure 1. The 2D density distribution of the errors between the 
redshifts determined by EaZY and those of[Ilbert et al.[ | [2009[ ) 
(COSMOS, Le Phare) and [Cirasuolo et al.| ( [2007[ l (UDS, Hyperz) 
for all narrow-band excess sources (full catalog). We find that ~ 
95% of our measurements are in agreement with the redshifts from 
the literature. The median error is measured as 0.038 without 
sigma-clipping. There are over-densities at 2 ~ 0.4, 0.8, and 1.5 
that conform to redshifts for our major emission-lines. We find 
outliers with errors up to ~ 2.5, but these are only 4% of our 
sample. In comparison to the wealth of spectroscopic data, we 
find a median photometric redshift error of A 2 /I -|- Zspec = 0.047 
in comparison to spectroscopic redshifts. 


color excess for bright narrow-band magnitudes (see fig. 3 
of Sobral et al.[[201^ I. 


2.3 Photo -2 Measurements 


We initially used the photo -2 measurements from |Ilbert[ 
et al. (2009) and Cirasuolo et al. (2007) for COSMOS and 


UDS, respectively, that were provided in the corresponding 
catalogs. The main problem of using those measurements 
is that in the UDS catalog, more than > 60% of sources 
are without photometric redshifts. This raises issues when 
selecting emitters via redshift selection. Although the color- 
color selection is effective alone in selecting H/I-|-[OHI] and 
[OH] emitters (see below and appendices for discussion), we 
wish to have emitters selected by two independent methods 
to act as a check-and-balances to robustly select emitters. 
Therefore, we measured the photometric redshifts for all of 
the UDS and COSMOS narrow-band excess sources, using 


EaZY (Brammer, van Dokkum & Coppi 20081 to ensure 


that (1) the majority of the UDS catalog had reliable pho¬ 
tometric redshifts, and (2) that both COSMOS and UDS 
had their redshifts determined by the same code, models, 
and assumptions. 

The filters we use for measuring photometric red¬ 
shifts for our UDS sources are UBVRizYJHK -I- Spitzer 
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Table 1. A list of the narrow-band filters used (Sobral et al. |2013t , along with the central wavelength (|J-m) and the FWHM (A) of each 
filter. Included is the expected redshift of each emission line within the range of the FWHM of the filter and the comoving volume that 
is observed. 


H^+[OIII]5Q07[OII]3727 


Filter 

■^obs 

(pm) 

FWHM 

(A) 

Z 

Volume 

(10^ Mpc^ deg“^) 

z 

Volume 

(10^ Mpc^ deg’ 

NB921 

0.9196 

132 

0.84 ± 0.01 

1.79 

1.47 ±0.02 

3.75 

NBJ 

1.211 

150 

1.42 ± 0.01 

3.11 

2.25 ± 0.02 

4.83 

NBH 

1.617 

211 

2.23 ±0.02 

5.05 

3.34 ±0.03 

6.53 

NBK 

2.121 

210 

3.24 ± 0.02 

4.87 

4.69 ±0.03 

5.68 






Figure 2. Log-scale photometric redshift distributions of the emitters selected in NB921{upper left), N'Bj{upper right), NBH(^ow^er left), 
and NBK(^otuer right). Each main peak is associated with a strong emission line, specifically Ho, H/3-|-[OHI], and [OH]. The dashed 
red line is the expected redshift of Ho emitters. The dashed blue lines are the expected redshifts for H/3-|-[OHI] emitters corresponding 
to H/34861, [OHI]4959, and [OHI]5007. It is clear then that differentiating these lines, even in narrow-band surveys, is quite difficult 
due to their close proximity to each other. Lastly, the dashed green line is the expected redshift of [OH] emitters. Highlighted in the 
corresponding, but lighter colors, are the photo- 2 : selection regions. The Ho selection region is from [Sobral et ah] ( |201^ . The photo- 2 : 
distributions are from our photo- 2 : calculations using EaZY for the COSMOS and UDS fields. 
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Table 2. The number of emitters selected based on color-color, photo- 2 :, spec- 2 :, and dual/multi-emitters. We also include the total 
number of emitters found. It should be noted that an emitter can be selected by more than just one selection, such that by tallying up 
the number of emitters per column will result in a number larger than the number in the column that shows the total number of emitters. 
We also show the number of sources selected only based on color-color and only based on photo-z. This highlights the importance of 
having more than one selection technique. For example, if we solely relied on color-color selection, then we would have a loss of 21% in 
our z = 1.42 H/3-|-[OIII] sample. Furthermore, we also include the number of emitters selected as dual/multi-emitters. These are sources 
that were detected in more than one narrow-band, resulting in two or more detected emission-lines complementing each other (e.g., Ha 
in NBH and [Oil] in NB921 for z ~ 1.47). 


Emission Line 

Band 

z 

Color-Color 

CC Only 

Photo -2 

Photo -2 Only 

Spec -2 

Dual/Multi-Emitter 

Total 

H/3-t[OIII] 

NB921 

0.84 

2005 

1000 

1262 

257 

213 

160 

2477 


NBJ 

1.42 

277 

41 

314 

78 

15 

23 

371 


NBH 

2.23 

208 

52 

212 

56 

3 

44 

271 


NBK 

3.24 

145 

11 

158 

24 

2 

0 

179 

[Oil] 

NB921 

1.47 

3152 

957 

2211 

16 

97 

213 

3285 


NBJ 

2.23 

115 

51 

85 

21 

0 

6 

137 


NBH 

3.30 

29 

18 

16 

5 

1 

0 

35 


NBK 

4.70 

18 

14 

4 

0 

0 

0 

18 


IRAC Chi -4-1- narrow-band (NB) -|- broad-band 
(BB) filters. For our COSMOS sources, we combine our 
UBVgRizJK+GALEX FUV & NUV+Spitzer IRAC Chi 


vey (Bradshaw et al. 2013 McLure et al. 20131. We also 


- 4 -l-NB-l-BB catalog with that of Ilbert et al. (20091. The 


benefit of this is that we include 12 Intermediate Subaru 
(e.g., IA427, IA464) bands and one Subaru NB711 band 
that were within the Ilbert et al. (20091 catalog. This re¬ 


sults in a total of 29 filters to constrain the measurements. 
Furthermore, the benefit of using the narrow-band filters in 
measuring the redshifts is that these filters specifically cap¬ 
ture emission-lines which are inherent in the SEDs. The Pe- 
gasel3 spectral library that is used in EaZY includes a pre¬ 
scription for emission-lines, which makes the measurements 
more accurate when including the narrow-band filters. Fig¬ 
ure 1^ shows the benefit of using narrow-band filters in the 
fitting process by the sharp peaks for which we expect the 
major emission lines to be located in terms of redshift. 

Figure shows the 2D density distribution of = 
(^EaZY — z)/0 + 0^ where z are the redshifts measured by 
Ilbert et al.| ( |2009 1 and Cirasuolo et al.| (20071, and the orig¬ 
inal photometric redshifts in the catalog. We find a median 
error of 0.038 (all sources without sigma-clipping). Figure 
also shows over-densities at z ~ 0.4, 0.8, and 1.5 which are 
expected as these are the most populated redshift slices in 
our sample, since they conform to major emission-lines. We 
find outliers up to cr^ ~ 2.5, but these only constitute a small 
fraction, such that ~ 95% of our sample are in agreement 
with [Ilbert et al.| ( |2009| ) and [Cirasuolo et al.j ( [2007[ ). 


2.4 Spectroscopic Redshifts 

We make use of the vast array of spectroscopic observations 
from the literature, which greatly enhances the reliability 
of our sample. In the COSMOS catalog, spectroscopic mea¬ 
surements are drawn from various studies as listed on the 
COSMOS website, as well as the zCOSMOS measurements 


from Lilly et al. (20071. The UDS catalog also includes mea¬ 


include FMOS measurements from Stott et al. (2013aI, 


DEIMOS & MOSFIRE measurements from Nayyeri et al.. 


in prep, PRIMUS measurements from Coil et al. (20111, and 


VIPERS measurements from Garilli et al. (20141. In total. 


surements from various publications that are highlighted 
on the UKIDSS UKIRT website, including the UDSz sur¬ 


we have 1269 emitters that have spectroscopic redshifts with 
661, 350, 177, and 81 emitters in NB921, NBJ, NBH, and 
NBK, respectively. This allows us to enhance the reliability 
of our sample and to test our photo-z and color-color selec¬ 
tions. In comparison to the spectroscopic redshifts, we have 
assessed the median errors of our photometric redshifts to 
be Az/1 -I- Zspec = 0.047 (without sigma-clipping). 


2.5 Selection of H/l-|-[OIII] and [Oil] Emitters 

The selection of potential H/3-|-[OIII] and [Oil] emitters is 
done by a combination of three different methods: (1) pho¬ 
tometric redshift; (2) color-color; and (3) spectroscopic red¬ 
shift. In this section, we will present the general selection 
criteria that we used to select our sample. For more detailed 
information about the specific selection cuts applied in each 
case, we refer the reader to appendix |A| 

With three different selection methods, conflicts can 
arise where one selection method provides a result that con¬ 
flicts with another method. To solve the issue, we prioritize 
the selection methods as such: 1) spectroscopic redshifts, 
2) photometric redshifts, and 3) color-color. If the emitter 
has a spectroscopic redshift, then it is selected based only 
on measurement. If it doesn’t have a spectroscopic redshift, 
then we select it based on its photometric redshift. Lastly, 
if the emitter has a photometric redshift but is not within 
the range of being selected as Hq (see photo-z selection cri¬ 
teria in Sobral et al. 2013| ), H/3-|-[OIII], or [Oil] or if the 
emitter does not have photo-z measurements, then it is se¬ 
lected based on the color-color criterion. In most cases, we 
find that emitters with photo-z within the redshift selection 
range are also found within the color-color selection area. In 
such cases, the emitters are selected based on both selection 
methods. 
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As shown in figure]^ the local peaks in the photometric 
redshift distributions are located around the expected red- 
shifts for H/3-|-[OIII] and [Oil] emitters. This signifies that 
we have many H/3-|-[OIII] and [Oil] emitters in our sample. 
We select our emitter candidates by defining a range in the 
distribution of Zphot that is centered on the expected red- 
shift of the emission line, which is aligned with the peaks 
in figure In most cases, the number of Ha emitters are 
the largest, followed by ]OIII] and [Oil]. In NB921 (fig¬ 
ure]^, H/3-|-[OIII] and [Oil] lines are the strongest, respec¬ 
tively, as these redshifts are near the peak of the cosmic 
star-formation history and also probe a much larger cosmic 
volume that Ha. Other populations of emission lines are 
found, such as Paschen series lines, Hel, and [SHI] (figure 

i- 

Color-color selections are also applied for each redshift. 
The selection criteria used for 2 ~ 1 — 3 are from |Sobral| 


et al. (20131 as these are in perfect agreement with our 


large spectroscopic sample, while the criteria used for the 
H/3-|-]OHI] and [OH] at 2 = 3.3 and [OH] at 2 = 4.5 are 
based on the dropout color-color selection known as the Ly¬ 
man break technique ( [Dickins^ 1998| [Stark et al. 20091. For 
some redshifts, we use more than one color-color selection to 
reduce the contaminations from lower and higher -2 sources. 
All color-color selection definitions are found in table iati 

Whenever available, we select sources based on their 
spectroscopic redshifts. Sources for which the spectroscopic 
redshift contradicts the photo -2 and color-color selection are 
removed from the sample. By including spectroscopically 
confirmed sources, we increase the size of our sample. We 
also note that we use all our spectroscopic redshifts to con¬ 
firm the robustness of our photo -2 and color-color selections. 

Also in the selection process is selecting dual/multi¬ 
emitters. For some sources that are selected as emitters with 
some emission line, there also exists another emission line 
in another band. For example, an emitter that is selected in 
NBJ as an H/3-|-[OHI] emitter (2 = 1.47) can also be selected 
in NB921 as an [OH] emitter if observed in that narrow- 
band. [Sobral et al.| ( |2012[ | used this same technique in a 
double-blind study to find [OH] emitters in NB921 by using 
selected Ha emitters in NBH as a proxy. The benefit to this 
technique is that it confirms emitters if they are selected in 
at least two bands; corresponding to getting a redshift out 
of two emission lines. In cases where we find dual/multi¬ 
emitters, we treat them as spectroscopic measurements as 
it is similar to having a spectroscopically confirmed emitter 
and include them in the sample. 

One major source of contamination that we may have 
is having selected a source as an [Oil] emitter when it is an 
H/3-|-]OHI] emitter and vice versa. Also, there are situations 
where a source is selected as one of the emitters of interest, 
but also falls into the color-color selection for another emis¬ 
sion line, or even Ha (as these color-color selections were 
used by Sobral et al.| ( |2009[ |2012[ 20131 to find such emit¬ 
ters). To overcome this degeneracy, we look at the photo -2 
distributions shown in figure as a probability distribution 
to assign the emission line based on the most probable line 
in the sample. In all cases, except for NB921, Ha is the most 
probable line. For NB921, H/3-|-[OHI] is the most probable. 

Another source of contamination is from misidentified 
lines. From the wealth of spectra that we have, the majority 
of misidentified lines are Ha. Further details on the types 


of misidentified lines can be found in appendix]^ Our total 
sample size is outlined in table 


3 LUMINOSITY FUNCTIONS AND 
EVOLUTION 

We use the traditional Vmax estimator to create our binned 
data. The binned data is defined such that: 


HL,) = 


AlogLfjj 


E 


log Lc 


-A log Lf, ,72. 


C{L,)V^ 


( 2 ) 


where Lj is the luminosity bin, A log Ltij is the bin-size, 
log Lc,j is the central luminosity of the bin, C{Li) b eing 


3.1 


the completeness of the source described in section 
and Vmax.i being the volume for which that source may be 
detected as described in section 


3.1 Line Completeness 


To assess the completeness, we follow the methodology pro¬ 


posed in Sobral et al. (2012 20131. We start with the full 


catalog that has all the emitters and includes sources which 
did not make our emitter selection. To measure the recovery 
fraction based on the emission line flux, we input a mock 
line flux starting at 10“^® erg s“^ cm“^ to all sources in the 
catalog. We then apply the excess selection criteria (E and 
EW cuts) used in Sobral et al. (20131 followed by our color- 
color selections. The recovery fraction is then defined as the 
number of sources recovered divided by the total number of 
sources in the catalog. This is then repeated after increasing 
the input mock line flux by small increments. The expected 
result is that at 10”'^® erg s“^ cm“^ the recovered fraction 
will be low and will increase as the input line flux increases. 
Figurej^shows the average completeness determined for our 
H/3-|-]OHI] and [OH] sources in the different narrow bands. 
The advantage of this technique is that we are not limiting 
our determination of the completeness to a certain model, 
but actually using the observed data itself to get the com¬ 
pleteness correction. 

There are some important points to be noted. First, 
these simulations are run separately per image. This is be¬ 
cause the depths of each image are not the same and thus 
can not be used together all at once to determine the com¬ 
pleteness correction. Secondly, we apply an uncertainty of 
20% of the completeness correction to the other uncertain¬ 
ties in quadrature in order to take into account the errors 
associated with this method of determining the corrections. 


3.2 Volumes Filter Profile Correction 

We calculate the volumes assuming a top-hat filter that has 
the same range as the FWHM of the actual filter. We re¬ 
port the probed comoving volume per square degree in ta¬ 
ble The volumes for each log-luminosity bins in the lu¬ 
minosity functions are reported in tables |B1| and |B2| for 
H/3-b]OHI] and [OH], respectively. 

Although a top-hat filter makes our calculations easy, it 
is not a true representation of the throughput of the filter, 
which requires us to apply a filter profile correction. There 
are two main effects that the filter profile correction takes 
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Figure 3. Average line completeness for the entire sample per 
band. Note that the completeness does vary between the COS¬ 
MOS and UDS fields, as well as from image-to-image within each 
field. This is because each image has a different depth. We com¬ 
pute the completeness based on each image to account for this 
discrepancy. 


into account. The first is the flux loss due to emitters that 
are close to the edge of the filter’s FWHM. Bright emitters 
at the wings will have a significant flux loss (close to 40%; 
depends on the filter) and would be considered as a faint 
source. This gives an overall bias in our sample population of 
faint sources and a lack of bright sources. The second effect 
is the volumes are corrected for the bright sources. Any faint 
source that is close to the wings of the filter will most likely 
not be in our sample, but bright sources will be detected 
as faint emitters. This then implies that our bright emitters 


cover a wider range of the filter, meaning a wider range in 
redshift, and, therefore, a larger volume. 


We use the method proposed in Sobral et al. (2009 


2012 20131. We correct for the filter by creating a mock 
sample of 10® fake emitters based on the luminosity function 
with the assumption of a top-hat filter. Random redshifts 
are assigned to each source and in a range covering the full 
filter profile, but not large enough that evolutionary effects 
of the luminosity function and cosmological structure biases 
the results. We assume a uniform redshift distribution. This 
mock sample will have the same distribution as the input 
top-hat luminosity function which we define as 4>Ty_. We 
then make a second mock sample with the same top-hat 
luminosity function but now apply the actual filter. This is 
done by: 


S]] Li„{z)T{z)dz 

i::nodz 


(3) 


where T{z) is the filter-response function in terms of redshift 
and Lin (z) is the luminosity of the emitter which is defined 
as Lin{z) = L ^(Zrand “ Zfliter), where Zrand is the randomly 
assigned redshift and zgiter is the matching redshift of the 
filter. This results in the loss of sources as some sources 
will have a redshift outside the range of the filter’s FWHM. 
The luminosity function from this population is defined as 
<(>Fiiter and is compared to <(>th in order to get the filter 
profile correction factor. The result shows that the bright 
sources are underestimated as expected, which changes the 
shape of the final luminosity function slightly. Specifically, 
it decreases the faint-end slope and increases the L* and <(>* 
slightly. This correction is applied to the LFs by dividing all 
the binned ^[L) data points by the correction factor. 


3.3 Luminosity Function Fitting 


There are several different functions that have been pro¬ 
posed in the literature to describe the observed luminosity 
function of the Universe (see | Johnston|20TT] for an in-depth 
review). We adopt the most-widely accepted Schechter func¬ 
tion to fit the observed luminosity function. The Schechter 
function is defined in its log-form as: 

( \ l + a 

logio L (4) 


where 4>* is the normalization of the luminosity function, L* 
is the characteristic luminosity, and a is the faint-end slope. 

We fit each luminosity function using the MPFIT fitting 
routine ( |Markwardt| |20d^ , which utilizes the Levenberg- 
Marquardt least-squares minimization technique to find the 
best-fit for a given function. For each fit, we take the best-fit 
values as our fitted parameters. For the la errors, we run a 
Monte Carlo simulation, similar to that of Ly et al. (20111. 
The simulation starts by selecting a random number that 
is drawn from a normal distribution that will perturb each 
data point, 4>(L), in the luminosity function within the Icr 
error bars. We also vary the bin size and center of the bin 
by perturbing the original bin size and bin center by a uni¬ 
form distribution. The fit is then run again and these steps 
are repeated for each iteration. A total of 10® iterations are 
done to get a probability distribution of the best-fit values 
from where the Icr error bars are then calculated. 
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Figure 4. Left: Presented are our H/3+[OIII] Luminosity Functions along with LFs from the literature. Included on the top horizontal 
axis is the logjg SFR that was derived via the [Osterbrock fc Ferland] ( |200^ calibration (see section [4.2| l. The darker data points are 
color coded to match the lighter LF fit {dashed lines). There is a clear evolution in the LFs up to z ~ 3. We find that our z = 1.42 and 
2.23 LFs are in reasonable agreement with the [OIII] grism spectroscopy study of [Colbert et al.] ( |2013^ at the bright-end, suggesting that 
we are selecting a reliable sample of [OIII] emitters. The major difference between our z = 0.84 and the z = 0.83 LF of jLy et al.| | |2007l l 
is probably due to sample size biases. Our sample is much larger, hence we are able to populate our brightest bins, causing a shift in L* 
to higher luminosities. Top Right: The evolution of (p^. from the H/3+[OIII] luminosity function. A strong, decreasing evolution is seen in 
(pi, from z = 0 to z ~ 3. This same evolution is seen by UV LF studies (see jOesch et al.|2010| for details). Bottom Right: The evolution 
of Li, from the H/3+[OIII] luminosity function. We see a strong, increasing evolution in L* up to z ~ 3. 


Our best-fit values are shown in table]^ We keep a fixed 
to a constant value of —1.6 and —1.3 for H/3-|-[OIII] and 
[Oil], respectively, as we are not able to fully constrain the 
faint-end. These values are drawn from looking at past work 
from previous narrow-band studies in order for our results to 
be comparable (e.g., H(?-|-[OIII], [Colbert et ar]|2013| [Oil], 


Bayliss et al. 20111. The drawback to this is that we are using 


low-z measurements of a as a proxy for the high-z universe, 
which can be an incorrect assumption. Low-z studies, such 
Colbert et al. ] ( [2013| ), [Ly et al^([2007[), and [Pirzkal et ^ 


(20131 for OIII ai id|Ly et al. (20071, Ciardullo et al. 


(20131, Takahashi et al.| ( [2007[ ); ? Bayliss et al. (20111, and 


Sobral et al.[ ( [2012[ ) for [Oil], have shown that the faint-end 
slope doesn’t evolve up to z ~ 1.5, while Ha surveys such 


Sobral et al. ( 2013[[^014p have shown no evolution in the 
faint-end slope up to z ~ 2.23. Furthermore, UV studies 
(e.g., Oesch et al. 2010 Smit et al. 20121, have shown no 
evolution up to z ~ 6 — 7. Based on these results, we keep a 


hxed as we constrain the bright-end rather than the faint- 
end. 


Lastly, the LF results are not corrected for dust extinc¬ 
tion, except when measuring the star-formation rate densi¬ 
ties. This is because many studies in the literature use very 
different extinction diagnostics, such that it becomes diffi¬ 
cult to compare various studies. Furthermore, our knowledge 
of the role of dust on emission-lines for the high-z universe, 
especially for the emission-lines of interest in our study, is 
still in development and requires future detailed investiga¬ 
tions. To simplify the use of our LFs by others in future 
studies, we present all the LF parameter results as uncor¬ 
rected for dust and AGN contribution. When discussing the 
SFRDs in section [^ we will include the results with and 
without dust and AGN corrections. 
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3.4 H/5+[OIII] Luminosity Function z ~ 0.8, 1.5, 
2.2, and 3.3 


We present here the results of the fitted Schechter func¬ 
tion to the H/3-|-[OIII] observed luminosity function out to 
2 ~ 3.3. This is the highest redshift determination of the 
H/3-|-[OIII] LFs currently to date and is the first time that 
the luminosity function has been constrained out to these 
redshifts. We present the results in figure]^ From z ~ 0.8 to 
3.3, we see a clear evolution in the shape of the LFs (figure 
1^ left). We also show on figure]^ the evolution of <()* and L* 
with a fixed to —1.6. It should be noted that there is a de¬ 
generacy between the fitted Schechter parameters, as shown 
in figure This needs to be borne in mind when interpret¬ 
ing the evolution of any single parameter, although figure 
indicates that our results are relatively robust. Based on our 
results, we find that 0* has been decreasing from 2 ~ 0.8 to 
~ 3.3. The opposite trend is seen in L* where it is increasing 
from 2 ~ 0.8 to ~ 3.3. 


Our results show a clear evolution in the luminosity 
function and are consistent with the same evolution seen in 
the results from the literature (Ly et al.|2007 Pirzkal et al. 


2013 Colbert et al.|2013| Sobral et al.|2015|. We report our 

-- ^ ^ ■ 


LF parameters in tablel^ We find that our 2 ~ 1.47 and 2.23 


LF agrees well with the LFs of Colbert et al. (20131 in the 


bright-end, but diverges at the faint-end. This matches with 
our discussion in the section below (see section [3.4.1[ ) where 
we predict the bright-end to be dominated by [OIII] emit¬ 


ters. The Colbert et al. (20131 study was part of the HST 
WISP program, covering 29 fields (0.036 deg^) in search of 
Ha, [OIII], and [Oil] emission line galaxies using WFC3 grim 
spectroscopy. Because this was a spectroscopic study and 
the fact that our LF matches (in the bright-end) with that 
of Colbert et al.j ( [2013 \ gives us confirmation that we are 
picking up the [OIII] emitters in our sample. The rise in the 
faint-end can then be attributed to the H/? emitters in our 
sample. 

In comparison to |Ly et alT] | |2007| ) , we see a clear devia¬ 
tion of the fits between our 2 ~ 0.84 and their 2 ~ 0.83. Al¬ 
though we do agree in terms of (()*, the main deviation in the 
LFs are in L* and a (was a = —1.44±0.09 in comparison to 
our —1.6 fixed faint-end slope). The study was based on deep 
optical imaging of the Subaru Deep Field (SDF) using the 
Suprime-Cam on the 8.2 m Subaru Telescope and was com¬ 
plemented with Subaru FOCAS and Keck DEIMOS spec¬ 
troscopy. The deviation could be due to biases from sample 
sizes such that our sample consists of more bright emitters 
to populate the bright-end, hence shifting L* higher. In com¬ 


parison to Sobral et al. (20151, for which this work was done 


in unison with, we find perfect agreement but our sample 
probes deeper by 0.1 dex. 

Figurej^shows the evolution of L* along with the results 
from other studies. There is a strong trend in which L* is 
increasing from 2 = 0 — 2.23 and then flattens. This trend is 


supported by Ly et al. (20071, Pirzkal et al. (20131, Colbert 


et al. (20131, and Sobral et al. (20151. Prior to this work, 
the 2 < 1 studies hinted to a rising trend in L*, which with 
our measurements and the 2 ~ 1.5 measurements of ISobraTl 
jet al.j ( |2015| ) has been confirmed up to 2 ~ 3. 

For the normalization of the LF, we see an evolution 
(figure]^ such that <()* drops as redshift increase up to 2 ~ 
3. This is consistent with the collection of UV LFs (i.e., 




41 42 43 44 45 

logioLfoiiii (ergs') 


Figure 5. Shown is the predicted [OIII] LFs from 2 : 0.8 to 2.2 

and compared to the 2 : ~ 0.71 zCOSMOS Type-2 AGN LF of 
[Bongiorno et aL](|2010||. The H/3 LFs are made by simply taking 
the Hq LFs of |SobraI et al.| | [^13[ | and assuming a fixed H^/Ha 
ratio to convert them. We find that our 2 : ~ 0.84 LF is [Olllj- 
dominated at the bright-end. Also, the level of AGN contribution 
is very little, except for log^Q L[oiii] > 43.0 erg s“^. 


Qesch et al.||2Q10 l, while our determination is based on a 


reliable H/3-|-[OIII] sample. Note that prior to this study, the 
H/3-|-[OIII] measurements in the literature paint the picture 
that the 0* evolution is flat up to 2 ~ 1. With the inclusion 
of our measurements, along with the 2 ~ 1.5 measurements 
of Sobral et al. (2015 I, we find that </>* strong decreases after 
2 > 1. 


3.f.l Predicting the [OIII] LF and AGN contribution 

The results highlighted above are for the H/3 and [OIII] emit¬ 
ters combined as one sample since we can not separate the 
two types of emitters based on photometry. For this, we 
need to conduct spectroscopic follow-ups to properly seg¬ 
regate the emitters. We attempted to separate the sample 


by using the Ha LF of Sobral et al. (20131. The advantage 
of using the LFs of [Sobral et al.j ( |2013| ) is that it is fully 
compatible since we are both using the same data set and 
methodology. We start by first removing Aua ~ 1.0 mag 
dust correction, then apply an H/3/Ha = 0.35 line ratio 
from Osterbrock & Ferland (20061 to get the observed H/3 


LF. We then applied Anp = 1.38 mag (based on Calzetti 
[et al.[ ( [2000[ ); assuming Aua = 1.0 mag) to the LFs and 
dust-corrected our H/3-|-[OHI] LFs using Ah/ 3 +[oiii] = 1-35 
(see section [4.3.2[ for details). The next step was subtracting 
our H/3-|-[OHI] LFs from the predicted H/3 LFs to get the 
predicted luminosity function for [OHI]5007 emitters. The 
results are shown in figure [^ 

We find that the [OIII] emitters in our sample com¬ 
pletely dominate the H/3-|-[OHI] LFs while towards the faint- 
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Figure 6. Left: Presented are the [Oil] Luminosity Functions along with those from the literature. The SFR calibration used to create 
the top horizontal axis is from |Kennicutt| ( |1998| l (see section |4.2[ | . The darker data points are color coded to match the lighter LF fit 
(dashed lines). We find that the evolution from the low- 2 : studies of [Gallego et al.] l |2002[ | and |Ciardullo et al.|l|2013[| to o ur 2 : = 4.7 LFs 
is quite strong and clear. We find that our 2 : = 1.47 LF is in agreement with the HiZELS [Oil] study of [Sobral et al.|(|2012|| an d the 
Subaru Deep Survey study ofE et al.|p007[ |. Our z = 2.23 is also in agreement with the CF-HiZELS study of |Sobral et al.| | |2015[ |. Top 
Right: The evolution in the normalization of the LF. We find that tj)* has been decreasing from z ~ 1.47 to z ~ 5. Bottom Right: The 
evolution of L*. We find a clear, strong evolution in L* all the way to z ~ 5. 


end the H/? emitters dominate. This is expected as the the- 
ratio is ~ 3 (for Z = 0.0004; |Oster- 
, which would segregate our sample 
such that the bright-end will be populated by [OIII] emitters 
and the faint-end with H/3 emitters. We also find an interest¬ 
ing feature where the normalization in the [OIII] LFs are the 
same, with the exemption of the z ~ 1.42 [OIII] LF which 
is slightly higher. This can imply that the relative contri¬ 
bution of H/3 is the same for all three LFs. We note that 
this is a qualitative assessment and subtracting a Schechter 
function by another Schechter function doesn’t result in the 
same functional form. 


oretical [OIII]/H/3 line 
brock fc Ferland||2006 


bin of Bongiorno et al. (20101 (log^Q L[oiii] ~ 43 erg s“^), 
but in disagreement for the lower luminosity bins. This im¬ 
plies that our brightest [OIII] emitters are primarily AGNs, 
but the fainter emitters are a combination of star-forming 
galaxies and AGNs, with the star-forming galaxies being the 
most dominant. Future spectroscopic follow-ups of our sam¬ 
ple would allow us to properly study the evolution of AGNs 
in the Universe. 


3.5 [Oil] Luminosity Function z = 1.47. 2.23, 3.3, 
and 4.7 


We also attempted to compare the z ~ 0.8 [OIII] LF 
with the zCOSMOS AGN Type-2 LF of [Bongiorno et al.| 
120101 to qualitatively assess the contribution of AGNs. The 
Type-2 AGN is the best candidate that could contaminate 
our sample as they have a continuum that is similar to nor¬ 
mal star-forming galaxies and they photo-ionize the same 
cold gas that is photo-ionized by hot massive stars. We find 
that we are in agreement only for the brightest luminosity 


We present here the results of the [Oil] luminosity function 
and the Schechter fit out to z ~ 5. The results are high¬ 
lighted in figure We see a clear evolution of the LF with 
redshift, with a large increase in the characteristic luminos¬ 
ity with redshift. The right-hand panels of figureshow the 
evolution of the fitted />* and L* parameters, and figure 
shows the degeneracy between the fitted values. 

Included on figure are data from the literature that 
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Bayliss et al.|2011 

Ciardullo et al. 

2013 

|Ly et al.||2007t [Takahashi et al.||2007t [Sobral et al. 

2012 

2015 

1 . We find that our 2 = 2.23 binned LF data is 


(20151. Because their sample size is ~ 4 times larger than 


our measurement, we have combined their LF data points 
with ours. The main effect is our measurement extends the 
combined LF 0.15 dex fainter. We note that the LFs ofISo-l 


bral et al. (20151 are directly compatible with our LFs as 


our study follows the same methodology. In fact, thejSobr^ 


et al. (20151 is specific to the NBJ determined LFs and the 


effects of cosmic variance while this study focuses on the 
evolution and extension of the LFs out to 2 ~ 5. 

Our 2 = 1.47 measurements are in perfect agreement 


with Sobral et al. (20121 and close to agreement with Ly 
|et al.| ( |2007[ ). We note that for |Ly et ah] ( |2007[ ) the faint- 
end was measured to be a = —0.78 ± 0.13 and —0.9 ± 0.2 


for Sobral et al. (20121 while we keep our faint-end slope 


hxed to —1.3 for all LFs. As seen in £gure[^ we find that 
we are not probing deep enough to fully s ee the turn in the 
LF for Z/[oii] < 41 erg s“^ as found by Ly et al. (20071 
and [Sobral et ah] ( |2012[ ). This is probably due to the fact 
that both studies used 2" apertures, while our study uses 
3” apertures (provided in the Sobral et al. (20131 catalog) 
to select 2 = 1.47 [Oil] emitters, which means that their 
studies are better at recovering faint emitters. 

We find a strong evolution in L*, as shown in £gure[^ 
for which a strong, rising trend is seen up to 2 ~ 5. We 


also see the same evolution in Ha studies (e.g., Sobral et al.| 
(20121), where L* is strongly increasing from 2 = 0 to 2. 


The same evolution of L* is seen in UV studies (Oesch et al. 


2010 ) up to the same redshift range. We notice some scatter 
for the low -2 studies (Ly et al. 12007 Takahashi et al.||2007[ 


[Bayliss et al.|2011[ [Sobral et al.|2012[ |Ciardullo et al.|2013p . 

This is primarily due to limitations in survey area and/or 
shallowness of the surveys. 

The evolution in />* is shown in figure along with 
measurements from the literature. Our measurements show 
a decreasing trend since 2 ~ 1.5 while for redshifts less than 
1.5 shows a flat evolution in (()*.The same evolution in (j)^, is 
also seen in Ha studies (e.g., Sobral et al. ( 2013| ) where after 
2 ~ 1 and up to 2 ~ 2, is shown to be decreasing. UV 
measurements ( Oesch et al.||20f0 ( also see a similar trend. 

Future wide surveys, such as Euclid and WFIRST, will 
be able to observe larger samples of emission-line galaxies 
such that our results can be used as predictions for such up¬ 
coming projects. By taking our LFs and integrating them 
to some flux/luminosity limit, these future surveys can es¬ 
timate the number of [Oil] emitters that can be detectable. 
Our LFs would then be quite useful as a tool to plan sur¬ 
veys studying [Oil] emitters out to 2 ~ 5. Furthermore, our 
luminosity functions can be used as a tool to gauge the level 
of low -2 interlopers in various studies, such as Lya studies 
at high- 2 . 


4 EVOLUTION OF THE STAR-FORMATION 
HISTORY OF THE UNIVERSE 

In this section, we present the star-formation history evo¬ 
lution of our ]OII] sample out to 2 ~ 5. We begin by mea¬ 
suring the level of AGN contamination and then present 


the calibrations used to get the star-format ion rate densities 
(SFRDs). We conclude with a discussion of the evolution of 
the SFRD based on [OH] emitters, the correction for dust, 
and our estimates of the stellar mass density evolution of 
the universe based on our SFRD fit. 


4.1 Contribution from AGN 


Active Galactic Nuclei (AGN) play an important role in the 
evolution of galaxies. Because AGN heat the cold gas that is 
photo-ionized by O & B-type stars in star-forming regions, 
the same emission-lines become present by both sources. It 
is then imperative that the SFRDs are properly corrected for 
AGN contamination to ensure that the sample is, by major¬ 
ity, a star-forming sample. Due to the low number-density 
of AGNs, it is difficult to use the current catalogs in the 
literature (e.g., Chandra-COSMOS) as a direct indicator on 
the level of AGN contribution/contamination to our sample. 
When comparing to Chandra-COSMOS ( [Elvis et al.||20(39'| ), 
we find 1, 0, 4, 0 for 2 ~ 0.84, 1.42, 2.23, and 3.24 for our 
H/3-|-]OIII] sample and 5, 2, 1, 0 for 2 ~ 1.47, 2.25, 3.34, and 
4.69 for our [Oil] sample. We also compared our catalogs 
to the XMM-COSMOS catalog ( Cappelluti et al.|2009 ( and 
found no matches. 

These matches themselves can not give us a complete 
indication of the level of our AGN contamination as they 
are X-ray flux-limited. We instead take advantage of the 
rest-frame 1.6p.m bump in the SEDs of star-forming galax¬ 
ies. This bump arises from the minimum opacity of H“ 
ions in the stellar atmospheres of cool stars. For AGNs, 
the bump is meshed in with various other emission (e.g., 
PAHs, silicate grains) resulting in a rising power-law SED 
after 1.6p.m in the rest-frame. We use the deep IRAC data 
in COSMOS and UDS and with the condition that redder 
colors are AGNs, signifying the rising SED after 1.6p.m, and 
anything with bluer colors are star-forming galaxies result¬ 
ing in the 1.6p.m bump. We measure the colors by taking 
the [3.6 - 4.5] > 0.1 (2 ~ 0.8), [4.5 - 5.6] > 0.1 (2 ~ 1.5), 
and [5.6 — 8.0] > 0.1 (2 ~ 2.2) for H/3-|-[OHI] emitters 
and [4.5 - 5.6] > 0.1 (2 ~ 1.47) and [5.6 - 8.0] > 0.1 
(2 ~ 2.23) for [OH] emitters. We find AGN contamina¬ 
tion for H,fl-b]OHI] is ~ 11.4%, ~ 18.5%, and ~ 28.8% 
for 2 ~ 0.8, 1.5, and 2.2, respectively. The amount of AGN 
contamination for [OH] is ~ 18.5% and ~ 19.4% for 2 ~ 1.47 
and 2.23, respectively. For 2 > 2.2 in [OH] emitters, we set 
the AGN contamination constant to that at 2 = 2.2 as this 
would require going beyond the last IRAC band. Note that 
these are upper limits for the level of AGN contamination 
such that our SFRDs are corrected for the highest contami¬ 
nation possible via the 1.6p.m bump technique. By compar¬ 
ing our SFRD measurements to other star-formation tracers, 
we can determine if the AGN correction was too high or not. 
But to reliably measure the correction will require follow-up 
spectroscopy of our sample to properly separate the AGN 
from the star-forming sample. We therefore apply our deter¬ 
mined AGN correction to the luminosity densities measured 
from the fully integrated LFs (decreases the SFRDs) and 
include 20% of the correction factor in quadrature with the 
luminosity density errors. 
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Figure 7. Shown is the interdependent evolution of L* and 0* for H/3+[OIII] {left) and [Oil] {right). The contours are color-coded 
to match the text color in the legend. The confidence levels are organized such that the darkest shade is the la level and the lightest 
shade being the 3cr level. There is a clear evolution that as redshift increases, drops and L* increases up to 2 ~ 3 and 2 ~ 5 for 
H/3-|-[OIII]and [Oil], respectively. 


Table 3. The Luminosity Function parameters and derived properties for all H/3-|-[OIII] and [Oil] emitters. Errors in (j!)* and L* are 
computed by running a Monte Carlo Simulation, displacing all the measurements of ^finai by la. The faint-end slope, q, was fixed as our 
data don’t go faint enough to constrain the faint-end properly. The luminosity density, p£, was calculated by taking the infinite integral 
of the LF. The SFRDs were calculated based on the |Kennicutt| | [l998[ | calibration ([Oil]) and jOsterbrock fc Ferland] l |200^ (H/3-|-[OIII]) 
with p*,comp being the completeness and filter profile corrected SFRD measurement. p*,corr is the completeness -|- filter profile -h dust 
corrected SFRD measurement. p*,AGN-corr is the completeness -h filter profile -h dust corrected AGN corrected SFRD measurement. 
We show every measurement as p*,comp is the most robust measurement, while the dust and dust -h AGN corrected measurements are 
based on assumptions regarding line ratios, dust extinction laws, and AGN selection methods. 
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4.2 Calibrations 


The star-format ion rate density is calculated via the lumi¬ 
nosity density from the LF at each redshift. The luminosity 
density is defined as: 

POO 

L= <E>(L)LdL = (?i*L*r(a-H 2) (5) 

Jo 


where jC is the luminosity density, is the normalization, 
L* is the characteristic luminosity, and a is the faint-end 
slope. Our determined luminosity densities are highlighted 
in table and consider the full range of luminosities. The 
star-formation rate density (SFRD) is then calculated by 
using the [Kennicutt (19981 diagnostics: 

psFR.oii = 1.4 X 10 "^^£[oii]A40yr ^Mpc ® (6) 


where a £[oii]/-CHa = 1.77 is assumed. We note that us¬ 
ing the [Oil] SFR calibration comes with several drawbacks, 
such as metallicity, reddening, and line ratio assumptions, 
but the [on] line is the brightest emission-line detectable 
at 2 > 1.5 where Ha falls in the infrared. We will present 
the effects of these drawbacks in a future study (Khostovan 
et al. in prep). Furthermore, we present the uncorrected for 
dust SFRD measurements to see, qualitatively, the evolution 
of the SFRD. This means that the results shown in figure]^ 
are lower-limits since any dust correction will just increase 
the SFRD measurements. 

We also use the derived relation of lOsterbrock fc Fer-I 


psFR,H^+oiii = 7.35 X 10 "‘^Tlp/s+oiii] A^oyr ^Mpc ® (7) 

to measure the H,a-b[OIII] SFRUj^ although this can not be 
taken as a purely star-forming indicator as there are sev¬ 
eral caveats behind it. We want to make this point specifi¬ 
cally clear; we do not use the Hfi-i- [OIII] SFRDs in fitting 
the star-formation and stellar mass assembly history of the 
Universe. We instead use it to compare the measurements 
to those of the [Oil] SFRDs and other tracers in the litera¬ 
ture to show if our sample is tracing a star-forming sample 
and whether or not if the H/3-|-[OIII] calibration is more of a 
“reliable” tracer of star-formation activity than previously 
thought. 


land ( |2006 


4.3 Star-Formation Rate Density Evolution of 
[Oil] Emitters 

4.3.1 Uncorrected-for-Dust 

Figure shows the evolution of the uncorrected-for-dust 
[Oil] SFR density for the first time and determined in a self- 
consistent way from 2 = 0 — 5. The evolution is clear and 
signihes that a peak that occurs at 2 ~ 2 — 3 and then there 
is a fall for higher redshifts. We include [Oil] measurements 
from the literature after uncorrecting them for dust and cor¬ 
recting the cosmology (the pre-2000 papers used non-ACDM 
cosmological parameters). These results plus ours can then 
be taken as a lower limit as any dust extinction correction 
would increase the SFR densities. We include a compilation 


^ We used the dust extinction curve of jCalzetti et aL] ( |2000| l for 
the H/3-|-[OIII] emitters and applied for all redshifts, such that 
^H/3+[Olll] = 1-35 mag (assuming A^a = 1-0 mag). 
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Figure 8. The uncorrected for dust SFRD evolution based only 
on [Oil] emission studies. We find that our 2 : = 1.47 and 2.25 LF 
continues the inverse power-law slope that is found from 2 : = 0 — 2 
in the majority of SFRD studies, and a continuous drop for 2 : > 2. 
This is the first time that [Oil] studies have gone beyond 2 : ~ 1.5 
in a reasonably, statistically constrained fashion. We find perfect 
agreement with |Sobral et ah] ( [2012} and are 0.15 dex off from [L^ 
|et al.| < |2007[l. Also, we find perfect agreement with the 2 = 2.25 
SFRD of |Sobral et al.| l |2015[ l . 


of SFR densities and LF parameters from various studies, 
using various diagnostics, all normalized to the same cosmol¬ 
ogy as that of this paper, same Kennicutt (19981 calibration 
with the same line ratio, and all with j4[oii] = 0 to make it 
easier for future studies to utilize. This compilation is found 
in appendix [c| We find that our 2 = 1.47 [Oil] SFRD mea¬ 
surement is in perfect agreement with [Sobral et ah (20121 
and Ly et al. (20071. We are also in perfect agreement with 
the measurement of Sobral et al. (20151 and our 2 = 2.23 
measurement. 

We note that Bayliss et al. ( 2012| made a measure¬ 


ment at 2 = 4.6 by observing the GOODS-S field using the 
NB2090 and Ks filters of the ESO HAWK-I instrument. The 
results of this study are restricted to a sample of only 3 gen¬ 
uine emitters. Although ground-breaking at the time, their 
SFRD estimate is severely limited by issues of sample size 
and cosmic variance. 


4.3.2 Dust & AGN Corrected SFRD 

To compare with other studies using different SFRD di¬ 
agnostics, we adopt a dust correction using the HiZELS 
2 = 1.47 measurem ent ofjHayashi et al.j ( |2013| ), Ana ~ 0.35 
mag. jHayashi et al.| ( |2013| ) studied Ha and [OH] emitters us¬ 
ing HiZELS data to conclude that the traditional Aua = 1.0 
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mag that has been used in the literature is overestimating 
the dust correction for [Oil] emitters at a = 1.47, such that 
these emitters are observed to have Ahc ~ 0.35 mag and 
are less dusty than previously thought. 

To test the dust extinction coefficient of|Hayashi et al.| 
poIs] ), we apply the traditional Ana ~ 1-0 mag to all four 
[Oil] SFRD measurements. We find that based on this dust 
correction, our measurements overestimate the Ha-based 


SFRD measurements of Sobral et al. (20131 and the radio- 


stacked measurements of Karim et al. 


(20111, which is im¬ 


pervious to dust extinction. The level of overestimation is 
such that our 2 : = 1.47 SFRD measurement and z = 2.23 
SFRD measurement was ~ 0.4 dex above the SFRD mea¬ 


surements of Sobral et al. (20131 and Karim et al. (20111. 


When using the Hayashi et al. (20131 dust extinction coef¬ 


ficient, we find that our SFRD measurements are perfectly 


matched with Sobral et al. (20131 and Karim et al. (20111, 
as seen in figure |9| 

We apply the Calzetti correction ( Calzetti et al.pOOO l 
with the [Hayashi et a'L] ( |2013| ) measurement of Aua ~ 0.35 
mag such that: 


7l[oii] ^ fc([OH]) 
Ana fc(Ha) 


( 8 ) 


where fc([OH]) = 5.86 and fc(Ha) = 3.31, resulting in 
A[oii] = 0.62 mag. We calibrate all the measurements to the 
same [Oil] SFR calibration of [Kennicutt (19981. All mea¬ 


surements hereinafter include AGN corrections as discussed 
in section KT\ 

Figure [^ shows our dust-corrected and AGN-corrected 
[OH] SFRD measurements. We also include a large compi¬ 
lation of studies from the literature which is a combination 


of the compilations of |Hopkins & Beacom 

(|2006|1, |Madau 

& Dickinson (20141, Gunawardhana et al. ( 

2013[), Ly et al. 


(20071, and our own compilation as a comparison (appendix 
0 “ We find that our measurements accurately reproduce 
the star-formation history of the universe up to 2 ~ 5. This 
is the first time that an [OH] study has ever accomplished 
such a measurement in a self-consistent manner. We find 
that the 2 = 1.47 and 2.23 perfectly agree with the HiZELS 
Ha measurements of Sobral et al. (20131 on figure]^ The 
AGN contamination in Sobral et al. (20131 assumed a sim¬ 
ple ~ 10%, which is backed by a detailed search of potential 
AGNs by [Garn et al. (20101. Based on the similarities be¬ 
tween our [OH] measurement and the independently AGN- 
corrected SFRD measurement of Sobral et al. (20131, we can 


conclude that the level of AGN contamination measured is 
reasonable and the methodology sound. 

Another key point is the stacked radio measurements of 
Karim et al. ( 2011| {pink squares on figure]^. The benefit 


of radio measurements are that they are impervious to dust, 
but have the downside of poor resolution and blending. We 
find that our 2 = 1.47 and 2 = 2.25 [OH] measurements are 
in agreement withjKarim et al!](|2011 (, such that thejHayashij 


et al. (20131 dust extinction coefficient does reliably correct 


our measurements to represent the dust-corrected SFRD of 
star-forming galaxies. 

We also find an interesting result when comparing the 
H/3-|-]OHI] SFRD measurements to our [OH] SFRDs and 
other measurements in the literature. As discussed above, 
the H/3-|-[OHI] calibration is considered in the literature as 
a “mixed” tracer of star-formation activity. Here we find that 


using the calibration of jOsterbrock fc Ferland {20061 with 
a Ah/3-i-[oiii] = 1-35 mag (based on the traditional Aua = 1 
mag), our H/3-|-[OHI] measurement for 2 = 0.84 matches 
perfectly with the Ha SFRD of Sobral et al. (20131 and the 
radio measurement of jKarim et al. (20111. The implications 
of this agreement shows that not only is the dust correction 
technique applied correctly, but also that the AGN correc¬ 
tion is accurate such that it is matching with the AGN- 


corrected Ha SFRD of Sobral et al. (20131. For our 2 = 1.43 


H/3-|-[OHI] SFRD, we find a perfect match with Karim et al. 
(|2011|), [Sobral et al. ( 2012[), [Sobral et al. (20131, and our 
[OH] SFRD. The 2 = 2.23 H/3-b[OHI] SFRD ma tches well 
with [Karim et al. ] ( |2011| , [Sobral et al.[ ( |2012| ), and our 
[OH] SFRD measurement. Lastly, we find perfect agree¬ 
ment between our [OH] SFRD and the H/3-l-[OHI] SFRD 
at 2 ~ 3.3. All these perfect agreements hint to the notion 
that the H/3-|-[OHI] SFR calibrations could in fact be more 
of a reliable tracer of star-formation activity than previously 
thought. Furthermore, this is also strong evidence to show 
that our H/3-|-[OHI] sample is dominated by star-forming 
galaxies and is a reliable sample. Also, our survey seems 
to be detecting H(?-|-[OHI] emitters that have more dust 
in comparison to [OH] emitters such that the traditional 
Ana = 1 rnag applies to the H/3-|-[OHI] sample and a lower 
dust correction applies to the [OH]emitters. This notion was 
proposed by Hayashi et al.[ ( |2013| for their ]OH] sample. 
Their conclusion was that dustier [Oil] emitters fall to lower 
luminosities that are below the detection limit, while the 
less dusty emitters, which will be apparently brighter, are 
detected. 

We fit the SFRD using our ]OH] SFRD measurements 


along with the [OH] measurements of Bayliss et al. (20111; 
Ciardullo et al.] ( [2013 1 and Sobral et al. (20121 to the 
parametrization of|Madau fc Dickinson (20141: 


logic P* 


(1 + 2 )*’ 


'! + [(! + 2 )/c]‘* 


M 0 yr ^ Mpc 


(9) 


where our fit results with a = 0.015 + 0.002, b — 2.26 + 0.20, 
c = 4.07 ± 0.51, and d = 8.39 ± 2.60. The fit is purely 
based on [OH] emitters, but we have also fitted for the cases 
of [OII]+Ha+Radio, [OH]+UV, and [OH]+Ha+Radlo+UV 
(see figure to show how our fit will vary based on the 
data that we use. Based on the [OH] fit, we see a drop 
at 2 > 3 that is slightly steeper than those determined by 
UV dropout studies (i.e., Bouwens et al.[[2011 2014 Oesch 


[et al.[[20T0 Schenker et al. [[2013 1. Despite this drop in our 

[OH] SFRD compared to the UV studies, we do find that 
the UV measurements are still within Icr. 

An important note to make though is that prior to this 
paper, there does not exist a study besides UV/Lya studies 
that have measured the SFRD up to 2 ~ 5 since 2 ~ 3. 
This is a crucial point since there has been no other study 
so far that could confirm the drop-out measurements, which 
are severely affected by dust extinction. Furthermore, this 
is the first time that the cosmic star-format ion history has 
been constrained based on a single tracer for larger volumes 
and up to 2 ~ 5. Our current measurements are the farthest 
that we can measure the [OH] SFRD due to the fact that the 
emission line would go past A-band and into the infrared. 
Future space-based narrow-band surveys, such as JWST and 
the Wide-field Imaging Surveyor for High-redshift {WISH), 
will be able to probe [OH] emitters up to 2 ~ 12, which 
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Figure 9. Our [Oil] dust & AGN corrected SFRD evolution with the [Oil] studies of |Bayliss et al.| < |2011[ |; [Ciardullo et al.|l|2013[|;|Sobral 
|et al.| ( [2013| l and [Sobral et alT] ( |2015| l, along with the results of this paper, that are used to fit the paranietrization of jMadau & Dickinson 
l |2014| |. The best fit is shown as the dashed line {dodger blue) and is only based on [Oil] measurements. We also include an extrapolation 
to higher-z {dashed-dotted turquoise line), as we don’t constrain this part of redshift space but can extrapolate based on our fit. The 
1-a region is highlighted in gold filled regions around the fit. The stacked radio study ofjKarim et al.|||201I|l and the Ha study of|Sobral| 
jet al.| | [20T^ are also shown as a comparison and are in agreement with our measurements. Our compilation of SFRD measurements (in 
gray) are a combination of our compilation and that of [Hopkins fc Beacom| ( |2006| l, [Madau Sc Dickinsonj | [^014| |, [Ly et al. j pOOTl l, and 
[Gunawardhana et al.||2013| |. We reproduce the SFRD evolution history of the universe based primarily on [Oil] studies with the peak 
of star-formation history occurring at z ~ 3. We also include the fits o f[Hopkins &: Beacom] | |2006[ l (IMF corrected to Salpeter) and that 
of [Madau fc Dickinson] | |2014[ |. We find that the [Hopkins fc Beacom] ([2006^ fit reaso nably matches our SFRD fit, while the [Madau fcj 
|Dickinson| | [2014[ | fits well until z > 2. This is mostly because the [Madau fc Dickinson] | [2014[ | fit is driven by the z > 5 UV measurements 
(which are not backed by spectroscopy), for which we do not include in our [Oil] fit. 


would allow us to compare and confirm the UV SFRD mea¬ 
surements at 2 > 5. 


SFRDs (e.g., Bouwens et al.pOll 2014 Oesch et al.||2010 


jSchenker et al.|2013|. 


We also compare our fit to those of |Hopkins fc Beacom| 
| |2006[ ) and [Madau fc Dickinson] ( |2014[ ) in figure 1^ For the 
z < 2 regime, we find that our [Oil] SFRD fit agrees well 
with all the other fits. For the z > 2 regime, we do see 
divergences based on the fit. In terms of the actual data 
points, we find that the Hopkins & Beacom (20061 is closest 
in agreement as it has a continuing SFRD up to a peak 
at z ~ 2.5 and a drop that continues through the high-z 
[Oil] measurements. The Madau & Dickinson (20141 is also 
in agreement for the high- 2 : measurements, but fails to match 
with the 2 ~ 2 — 3 peak. This is mostly due to the fact that 
their measurements are driven by the 2 > 5 UV dropout 


As with all SFR measurements, there are systematic 
uncertainties that must be taken into account. In the case 
of [on] emitters, our main systematic uncertainties come 
from metallicity and dust extinction. To study the metallici- 
ties and its effects on the star-formation rate calibration, we 
will need to conduct follow-up spectroscopy. Furthermore, 
studying the metallicity of our sample will give us also an 
understanding of the dynamics (inflow/outflow) that can af¬ 
fect star-formation activity. We also plan to study in a fu¬ 
ture paper the dust extinction properties of our sample and 
how it relates to and affects the star-formation activity of 
galaxies in our sample (Khostovan et ah, in prep). 
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4.4 Evolution of the Stellar Mass Density 


We use the [Oil] SFRD results presented in this paper to 
provide an estimate of the stellar mass density (SMD) evo¬ 
lution by doing a time-integral of equation]^ The SMD evo¬ 
lution gives us an understanding of how the universe has 
assembled its mass throughout cosmic time. This estimate 
is quantitatively sensitive to the choice of the IMF in terms 
of the normalization of the SFRD and SMD evolution, but 
it does not qualitatively affect the final results. 

Our estimate assumes a Salpeter IMF, which has been 
used throughout this entire paper, and a recycling fraction 
of i? = 0.27. For a review of the derivation of this factor, we 
refer the reader to the recent review of IMadau fc DickinsonI 


(20141. We calculate the SMD by: 

I 


p 4 z) = {1-R) 


p*iz) 


Ho{i -b z)^nM{i + z)^ + nA 


measurements from Arnouts et al. (20071, Eisner, Feulner 

& Hopp (20081, Gallazzi et al. 

(20081, Perez-Gonzalez et al. 

(|2008|), [Kajisawa et al.|(|2009l 

Li & White] (|2009D, 

March- 

esini et al. (20091, Yabe et al. 

20091, Pozzetti et al. 

(20101, 


o 

ei, 


6X) 

o 


dz (10) 


where p*{z) is the SFRD fit using the parametrization 
defined in equation and R is the recycling fraction, or the 
fraction of stars that is returned back into the ISM and IGM. 
Because the equation above is using 2 : = oo as a reference 
for which we do not know the SMD for, we instead constrain 
the integral such that at 2 = 0 the SMD will be logj^Q p* ~ 
8.6 M 0 Mpc“®, in agreement with measurements made at 
that redshift. 

Our results are shown in figure |10[ We find an evolu¬ 
tion where the stellar mass assembly rapidly increases from 
10®'^ to 10® Mq Mpc“® from 2 ~ 5 to 2, a time frame of 
only 2 Gyr. The evolution then tapers and flattens out by 
2 = 0 which is related to the decrease in the SFRD that 
we have observed since 2 ~ 2. This is also the same conclu¬ 
sion found by observational studies of the SMD. We include 


2012 


2013 
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Caputi et al. (20111, Gonzalez et al. (|2011|, Bielby et al. 


Moustakas et al. 


(20131, and Muzzin et al. (20131 in 


the same evolution seen by these studies. We have found 
that the [Oil] based SFRD and SMD accurately reproduce 
the evolution of mass assembly in the Universe. This match 
can also be seen as yet another verification that our sample 
of [Oil] emitters are primarily star-forming galaxies as the 
conclusions from our SMD estimate are the same seen in the 
literature. 


5 CONCLUSIONS 

We have presented the largest sample of H/3-|-[OIII] and 
[Oil] emitters between 2 ~ 0.8 — 5 that have been selected 
based on a robust and self-consistent technique, backed up 
by a wide array of spectroscopic emitters. We have used the 
HiZELS UKIRT and Subaru narrow-band catalogs, along 
with multi-wavelength data from the COSMOS and UDS 
fields, to create a clean and well-defined sample of star¬ 
forming galaxies. The main results of this paper are as fol¬ 
lows: 

(i) We have robustly selected a total of 2477, 371 , 270, 
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Figure 10. The evolution of the stellar mass density of the Uni¬ 
verse based on the integrated [Oil] SFRD. Overlaid are the SMD 
measurements from the literature that were compiled in the recent 
review of |Madau fc Dickinson| | |2014[ |. We find that our integra¬ 
tion of the purely [Oil] determined SFRD reasonably traces the 
stellar mass assembly of the Universe. 


179 H/3-l-[OIII] emitters at 2 = 0.84, 1.42, 2.23, and 3.24 and 
3285, 137, 35, 18 [Oil] emitters at 2 = 1.47, 2.25, 3.34, and 
4.69 in the combined COSMOS and UDS fields. These are 
the largest samples of H/3-|-[OIII] and [Oil] emitters to have 
been detected in this redshift range. 

(ii) We have extended the luminosity function in the 
literature to higher- 2 , as well as refined the lower -2 measure¬ 
ments for both types of emitters. For the H/l-|-[OIII] emit¬ 
ters, we find that the bright-end of our 2 = 1.42 and 
2 = 2.23 LFs are in agreement with the grism spectroscopy- 
based luminosity functions of Colbert et al. (20131; hence, 
this increases the reliability of our sample being dominantly 
[OIII] emitters in the bright-end. We also find from our pre¬ 
dictions of the [OIII] LFs that our sample is dominated by 
[OIII] emitters at the bright-end. The faint-end is dominated 
by H/3 emitters. We also find that the normalization of the 
[OIII] LFs are the same such that the relative contribution 
of H/3 emitters is the same between 2 ~ 0.8 — 2.2. (iii) 
The evolution of L* and for H/3-|-[OIII] is found to have 
a strong increasing/decreasing evolution, respectively, up 
to 2 ~ 3. For our [Oil] sample, we find that L* increases 
strongly up to 2 ~ 5 and 0* is strongly dropping up to the 
same redshift. 

(iv) We have discussed that our luminosity functions 
are reliable to be used in making predictions of the num¬ 
ber of emitters to be detected by future wide-surveys, such 
as Euclid and WFIRST. Furthermore, our luminosity func¬ 
tions can also determine the number of low -2 interlopers in 
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LyQ studies, such that the level of contamination by low -2 
sources can be reduced in such studies. 

(v) The SFRD has been constrained using [Oil] mea¬ 
surements up to 2 ~ 5 for the first time. We find that the 
peak of the cosmic SFRD is located around 2 ~ 3 and is in 
agreement with our large compilation of UV, IR, radio, and 
nebular emission studies. We find that for 2 > 2, our SFRD 
fit drops slightly faster in comparison to the UV dropout 
studies in this redshift regime. However, we find that the 
UV measurements are within the 1-a error bar range of our 
SFRD fit. Future space-based narrow-band surveys, such 
as JWST and WISH, will be able to extend the range of 
[Oil] detection out to 2 ~ 12 so that we can compare and 
confirm or invalidate the UV dropout measurements. 

(vi) We also find that the H/3-|-[OIII] SFRD measure¬ 
ments are nicely in line with our [Oil] sample and other star- 
formation tracers. This then brings to question of whether 
the H/3-|-[OIII] calibration is more “reliable” as a tracer of 
star-formation than previously thought. With our large sam¬ 
ple of these emitters, we will have the ability to explore this 
issue in detail. 

(vii) By integrating the SFRD, we have made estimates 
of the stellar mass density evolution and find that it steeply 
rose up to 2 ~ 2 and flattened out up to the present-day. 
This is also confirmed by the wealth of measurements in the 
literature. 

The results in the paper have implications in the evo¬ 
lution of galaxies and the star-formation activity occurring 
in said galaxies. Despite the robustness of our sample, there 
is still room for improvement. Our measurements have done 
well to constrain the bright-end, while keeping the faint-end 
fixed based on measurements from the literature. We will 
require deeper narrow-band and broad-band measurements 
in order to constrain the faint-end slope of the LF. Spectro¬ 
scopic follow-up will also be necessary to accurately measure 
the extent of AGN contamination in our sample. Although, 
our color-color selections have shown (see figure A11 that 
they are quite reliable due to the large set of spectroscopic 
measurements confirming this reliability. That being said, 
spectroscopic measurements of our sample will help in sepa¬ 
rating the H/3 and [OIII] samples to measure separate lumi¬ 
nosity functions. Lastly, future narrow-band surveys, such as 
the proposed VF/Si/telescope, will be able to extend the red- 
shift window of H/I-|-[OIII] and ]OII] studies up to 2 ~ 12, 
which can be used to confirm the UV dropout studies at 
higher- 2 . Despite all these improvements and potential fu¬ 
ture progresses, our sample has reliably (given all the lim¬ 
itations) and robustly traced the evolution of star-forming 
activity in the universe. 
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Figure Al. The i — K versus B — R color-color distribution of 
NB921 emitters using the color-color selection of |Sobral et ah] 
l |2013| |. Overall, the level of completeness is > 90% as the large 
majority of spectroscopically confirmed emitters are within the 
selection area; hence, we confirm the the BRiK selection of |SobraT| 
|et al.| ( |2009| l is very efficient in selecting Ha, H/3+[OIII], and 
[Oil] samples from narrow-band surveys. 


APPENDIX A: SELECTION TECHNIQUE 

Here, we present, in detail, the selection of emitters that 
made it in to our sample. We also present the exact color- 
color selections that were applied in our work in table [XT] 


AO.l UP+IOIII] Emitters at z ^ 0.8 


H/3-|-[OIII] sources in the NB921 data at z ~ 0.8 are selected 
by their photometric redshifts within the range of 0.75 < 
Zphot < 0.95. The color-color selection criterion reduces the 
number of contaminants by separating the H/3-|-[OIII] emit¬ 
ters from lower -2 Ha and higher -2 [OH] emitters. This is 
done by using the BRiK selection (Sobral et al. 20091 in 
hgure I 


IT] Spectroscopic redshifts were used to assess the 
robustness and effectiveness of the selection criteria. 213 
sources were spectroscopically confirmed. 169 were selected 
by the color-color selection. From the 213 sources (for which 
all were selected by their photo- 2 ), only 11 were confirmed 
H;g4861, 76 were [OHI]4959, and 126 were [OHI]5007 emit¬ 


ters. Removed from the sample where 9 low -2 and 13 high -2 
spectroscopically confirmed emitters. The lower -2 emitters 
were primarily Ha and [Nil]. The higher -2 emitters were pri¬ 
marily [OH] emitters. All these misidentified emitters were 
removed from the sample. Based on the spectroscopic data, 
we find that ~ 91% of all spectroscopic measurements for 
photo -2 selected objects were either H/3 or [OIII] emitters 
and that the BRiK color-color selection does select ~ 91% 
of all the H/3-|-[OHI] spectroscopically confirmed emitters. In 
total, we have 2477 2 = 0.84 H/3-|-[OHI] emitters in our sam¬ 
ple. Based on the spectroscopic data, we find that the color- 


color selection effectively selects H/3-|-[OHI] emitters with 
a completeness of ~ 91%, making our sample not just the 
largest, but the most complete sample of H/3-|-[OHI] emitters 
to date. 


AO.2 HP+[OIII] & [Oil] Emitters at z ^ 1.5 

[OH] emitters in NB921 are selected with photometric red¬ 
shifts between 1.2 < ZpUot < 1-7. We use the BRiK color- 
color selection and include any sources with spectroscopic 
redshifts. Included in our sample are 97 spectroscopically 
confirmed sources, with 90 of them being color-color selected 
as well. Removed from the sample were 48 low -2 spectroscop¬ 
ically confirmed emitters, for which the majority were [OIII]. 
We also removed a few [Hel] emitters at 2 ~ 1.28 and Ha, 
and [Nil] emitters (2 ~ 0.4). The contamination from high -2 
emitters was significantly less (8 emitters) as there are no 
major emission lines beyond [OH]. The majority of this con¬ 
tamination came from [Mgll] emitters at 2 = 2.25. The issue 
of contamination arises here as we may “naively” state that 
our level of contamination is ~ 33%, but it is noted that 
spectroscopic measurements to date have an inherent bias 
to the low -2 regime, such that there are more low -2 than 
high -2 measurements. This makes accurately measuring the 
level of contamination difficult. We note though that our 
color-color selection did select ~ 93% of the spectroscop¬ 
ically confirmed emitters. In total, we have selected 3285 
2 = 1.47 [OH] emitters. This is by far the largest sample of 
[OH] emitters at 2 ~ 1.5 to date and, based on the spectro¬ 
scopically confirmed sources, is ~ 93% complete. 

H/l-|-[OHI] emitters in NBJ are selected based on a pho¬ 
tometric redshift range of 1.20 < Zphot < 1.70. The color- 
color selection criteria consists of a BzK and izK selection, 
as shown in figure |A2[ We use the BzK selection to get 
our initial sample of emitters and remove the lower -2 con¬ 
taminants, which are mostly Ha emitters. To remove the 
higher -2 contaminants ([OH]), we use the izK selection. 
There were also 15 spectroscopically confirmed sources, all of 
which were within our color-color selection. Of these 15 emit¬ 
ters, 4 were H/3, 5 [OHI]4958, and 6 [OHI]5007 emitters. We 
removed 5 emitters that were spectroscopically confirmed. 
These emitters were primarily Ha and [Nil], all of which 
were removed from the sample. Our final sample consists of 
371 H/3-|-[OHI] emitters at 2 = 1.47 that were selected. 

A0.3 HP+[OIII] & ]OII] Emitters at z ^ 2.2 

[OH] emitters in NBJ are selected if their photometric red¬ 
shifts are between 1.7 < Zphot < 2.8. We apply the BzK 
color-color selection to remove the lower -2 contaminants, 
which are primarily Ha emitters, as shown in hgure [A^ We 
then use the izK color-color selection to separate the sam¬ 
ple from H/3-|-[OHI] emitters. There were no spectroscop¬ 
ically conhrmed sources included in the sample. Removed 
from the sample were 3 contaminants, which were all low -2 
emitters (1 [Nil], 1 [OIII], and 1 [Hel]). In total, there are 
137 [OH] emitters selected at 2 = 2.25. 

H/3-l-[OHI] emitters at 2 ~ 2.2 in NBH are selected with 
photometric redshifts between 1.7 < Zphot < 2.8 along with 
BzK and izK color-color selections. We apply the BzK 
color-color selection to remove emitters with 2 < 1.5 (pri¬ 
marily the 2 = 1.47 Ha emitters) and then apply the izK 
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Table Al. Definitions of the Color-Color Selection 


Filter 

Color-Color 

Emitter 

Redshift 

Selection Criteria 

NB921 

BRiK 

H/3+[OIII] 

[Oil] 

0.84 

1.47 

-0.3 <{B -R) < 0.08 h {i- K) < 2.04(B - R) + 0.81 

0.08 <{B- R) <l 1.6(5 - R) <{i- K) < 2.04(5 - R) + 0.81 

1 < (5 - 5) < 1.24 & 1.6(5 - R) <(i- K) < 3.21 

1.24 < (B- R) k 2.01 <{i- K) < 3.21 
(B-R) < -0.3 

-0.3 <{B- R) < 1.17 k {i- K)> 2.04(5 - 5) -I- 0.81 
(5 - 5) > 1.17 k (i-K)> 3.21 

NBJ 

BzK 

H/3-|-[OIII] & [Oil] 


{B -z) < 0.4 

0.4 <(B -z) < 2.41 k (z - K) > (B - z) - 0.4 

2.41 <{B - z) k (z - K) > 2.0 


izK 

H/3+[OIII] 

[Oil] 

1.42 

2.25 

(2 — K) < 5(2 — 2 ) — 0.4 
(2 — K) > 5(2 — 2 ) — 0.4 

NBH 

BzK 

H/3+[OIII] 

2.23 

(B -z) < 0.4 

0.4 < (5 - z) < 2.41 k {z - K) > {B - z) - 0.4 

2.41 < (B - z) k {z - K) > 2.0 


izK 

H/3+[OIII] 

2.23 

(2 — K) > 5(i — z) — 0.4 


UVz 

[Oil] 

3.34 

{U -V)> 1.2 k(U -V)> 0.5(y - 2 ) + 1.2 k(V -z) < 1.6 

NBK 

UVz 

H/3+[OIII] 

3.24 

{U -V)> 1.2 k{U -V)> 0.5(y - 2 ) + 1.2 k{V -z)< 1.6 


Viz 

[Oil] 

4.69 

{V -i)> 1.2 k{V -i)> 0.89(1 - 2 ) -I- 1.2 k (i - z) < 1.3 
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Figure A2. Color-color magnitude distributions for all NBJ emitters. Left: The z — K versus B — z selection used to separate the 
H/3-|-[OIII] emitters at 2 ~ 1.47 and [Oil] emitters at 2 ~ 2.23 from the low -2 emitters that are primarily Ho. Right: The z — K versus 
2 — 2 selection is used to select [Oil] and H/3-|-[OIII] emitters. For both selections, the spectroscopically confirmed emitters lie within 
the selection region adding to completeness of the sample. Both color-color selections nicely distinguish between the two samples. We 
note that for the izK selection, about ~ 15% of selected [Oil] emitters are within the selection region of H/3-|-[OIII] emitters. These 
are photo -2 selected and shows that relying on purely the color-color selection would show a ~ 15% drop in the completeness of the 
[Oil] sample, and a ~ 15% increase in the contamination of the H/3-|-[OIII] sample. 
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Figure A3. Color-color magnitude distributions of all NBH emitters. Left: The z — K versus i — z color-color selection used to separate 
the H/3-|-[OIII] 2 : 2.23 emitters from the Ha 2 : ~ 1.47 emitters (selected based on the methodology of |Sobral et al.||2013[ |. We have 

highlighted the spectroscopically confirmed emitters. The color-color selection shows a clear separation between H/3 and Ha measurements 
based on spectroscopic measurements alone. We find that 21% of our photo- 2 : selected H/3-h[OHI] emitters are within the selection area 
for Ha emitter, showing that relying on just color-color selection would result in the loss of ~ 56 emitters from the sample and an increase 
in contamination of the Ha sample. Right: The U — V versus V — z color-color selection that is based on the Lyman break drop-out 
technique and is used to find 2 ~ 3.3 [OH] emitters. We also include the H/3-|-[OHI] sample and show that the vast majority of these 
emitters are outside the [OH] color-color selection region. 
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Figure A4. Color-color magnitude distributions of all NBK emitters. Left: shows the U — V versus V — z color-color selection used 
to find 2 ~ 3.3 H/3-|-[OHI] emitters based on the Lyman break drop-out technique. We also include the Ha emitters (2 ~ 2.23) and 
find that this selection reliably separates the two samples. Right: shows the V — i versus i — z color-color selection used to find 2 ~ 4.7 
[OH] emitters based also on the Lyman break drop-out technique. We include the H/3-|-[OHI] sample to show that this selection criteria 
is reliable to separate the two samples. We find no H/3-|-[OHI] emitter falling within the [OH] color-color region. 
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color-color selection (shown in figure A31 to separate our 
H/3-|-[OIII] sample from the z = 1.47 Ha emitters. We find 
that 21% of our sample is outside the color-color region in 
figure but is still selected via the photo -2 selection. This 
would raise concerns about the reliability of the color-color 
selection, but we must point out that it is reliable in terms of 
the separation between the H/3-|-[OIII] and Ha spectroscop¬ 
ically confirmed emitters. This also raises the point that we 
are more concerned with consistency in the use of selection 
techniques to reduce the effects of assumptions in our sam¬ 
ple. Furthermore, this also shows that we can not rely on the 
color-color selection technique alone to select emitters since 
this would result in a 21% drop in H/3-|-[OHI] emitters and a 
21% increase in the contamination of the Ha sample in|So-| 


bral et al. (20131. By including the photo -2 in our selection 


methodology, we reduce the level of contaminants and also 
increase the reliability of the sample. We also included three 
spectroscopically confirmed sources in the sample that were 
selected by their color-color selection. Of these 3, we have 
two [OIII]4959 and one [OHIJ5007 emitter. Removed from 
the sample were 3 low -2 spectroscopically confirmed emit¬ 
ters (Ha at 2 ~ 1.47). There were no high -2 spectroscopic 
measurements that contaminated the sample. A total of 271 
H/3-|-[OIII] emitters at 2 = 2.23 were selected. 


AO.4 Hl3+[Oin] & [Oil] Emitters at z ^ 3.3 

[OH] emitters at 2 = 3.3 in NBH are selected if 2.8 < Zphot < 
4. We also select sources if they satisfy the UVz color-color 
criteria. This separates the lower -2 contaminants from our 
sample. We include our H/3-|-[OHI] sample in figure [A3| to 
show the separation between the [OH] and H/3-l-[OHI] sam¬ 
ple and find that the UVz selection is reliable in selecting 
our [OH] sample. Furthermore, as there are no major emis¬ 
sion lines that are detected at higher- 2 , there is no need 
to include another color-color criteria to account for this 
contamination. As shown in fig. the number of sources 
greatly drops at higher -2 for the redshift of interest, making 
the number of contaminants very small. There was only one 
spectroscopically confirmed source that was included in the 
sample and only found in COSMOS. The UVz color of this 
emitter places it well within the selection area, adding to the 
reliability of our color-color selection. No spectroscopic con¬ 
firmed emitters were misidentified in the selection. In total, 
there are 35 of [Oil] emitters at 2 = 3.3 selected. 

Our H/3-|-]OIII] emitters at 2 = 3.3 in NBK are selected 
if their photometric redshifts lie between 2.8 < Zp^ot < 4. 
We use a UV z selection criteria, as shown in figure [Ad] to se¬ 
lect H/3-|-[OIII] emitters based on the Lyman break dropout 
technique. Sources with no detection bluer than the [/-band 
and detection in the V and 2 bands greater than the 5a mag¬ 
nitude detection limits were included. We also include all Ha 
emitters with photo -2 around 2 = 2.23 in figure [AT] to show 
the separation between them and our H/3-|-[OIII] selected 
emitters. One spectroscopically confirmed source was also 
included in the sample from UDS. We find that this spec¬ 
troscopically confirmed emitter is well within the color-color 
region. No spectroscopic confirmed emitters were misidenti¬ 
fied in the selection. A total of 179 H/J-|-[OHI] emitters at 
2 = 3.3 were selected. 


AO.5 [Oil] emitters at 2 ~ 4.7 


We select our [OH] emitters in NBK if they have a photo¬ 
metric redshift between 4.0 < Zphot < 6.0. Emitters are also 
selected by using the 1^/2 criteria of Stark et al. (20091 for V- 


band dropouts. The color-color selection is shown in figure 
A4 and includes our H/3-|-[OHI] sample. We find that there 
are no H/?-|-[OHI] emitters within the [OH] selection area, 
which adds to the reliability of our sample. We find only 
one emitter was selected by its photo -2 within the selection 
range and it is well within the color-color selection region. 
The majority of emitters were selected by the color-color se¬ 
lection. All sources selected were also under the condition 
that anything bluer than F-band must have no detection. 


A study of [OH] emitters by Bayliss et al. (20121 also used 
a similar technique using the BVz criteria of [Stark et al.| 
(20091 for their 2 ~ 4.6 sample. It must be noted though 


that the results of this study should not be taken as reliable 
due to the fact that they were limited to a sample size of 
only 3 [OH] emitters (about 3 times smaller than our sam¬ 
ple) and the volume probed by this study is a factor of 100 
times smaller than our study making it severely susceptible 
to cosmic variance. Our sample is statistically larger and ro¬ 


bust in comparison to Bayliss et al. (20121. There were no 


spectroscopically confirmed sources in our sample and no 
spectroscopic confirmed emitters were misidentified in the 
selection, thus, giving a total of 18 [Oil] emitters at 2 = 4.5 
that were selected. 


AO.6 Notes on Contamination 

We advise the reader that the measurements of contami¬ 
nation are not strictly reliable for the H[?-|-[OHI] 2 ~ 1.5 
emitters due to the bias in the spectroscopic redshift dis¬ 
tribution. Two issues arise are: (1) lack of spectroscopic 
measurements at higher redshifts to properly quantify the 
level of contamination, and (2) the inherent bias of spectro¬ 
scopic measurements to the lower -2 regime. The first point 
really just requires more spectroscopic measurements to in¬ 
crease the population of spectroscopically confirmed sources. 
The second point has to do with the distribution of spec -2 
measurements. There exists more spec -2 measurements for 
2 < 1, which results in a skewed histogram that favors the 
lower -2 regime. When measuring the level of contamination, 
there are more low -2 measurements than spectroscopically 
confirmed measurements (for example, 2 = 1.47 [OH]) which 
causes a “naive” and biased estimation of the level of con¬ 
tamination. Instead, we considered where the spectroscopi¬ 
cally confirmed H/J-|-[OHI] and [Oil] measurements were on 
the color-color diagrams as a way to assess the reliability of 
our selection technique. 

Based on the points described above, we can mea¬ 
sure the level of contamination for the H/3-|-[OHI] 2 ~ 
0.84 sample (~ 10%). To robustly measure the contami¬ 
nation for the higher -2 samples, we will need to conduct 
spectroscopic follow-up which is currently underway with 
Keck/MOSFIRE and ESO/VLT (Khostovan et ah, in prep). 
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Table Bl. H/3+[OIII] Luminosity Function. <I>oijs shows the ob¬ 
served LF data points per bin; simply, it is the log of the number 
of emitters divided by the volume. 'Ffinal is the completeness and 
filter profile corrected luminosity data points per bin. The errors 
here are Poissonian but with 20% of the corrections added in 
quadrature. 


Table B2. [Oil] Luminosity Function. 3>oba shows the observed 
LF data points per bin; simply, it is the log of the number of 
emitters divided by the volume. <I>flnal is the completeness and 
filter profile corrected luminosity data points per bin. The errors 
here are Poissonian but with 20% of the corrections added in 
quadrature. 


loglO ^H/3-|-[OIII] 
(erg s“^) 

# 

^obs 

(Mpc“® dlogio L) 

^final 

(Mpc“3 d logic L) 

Volume 
(10® Mpc^) 

z = 0.84 

41.10 ±0.10 

703 

-1.97 

-1.82 ± 0.02 

3.25 

41.30 ±0.10 

465 

-2.15 

-2.04 ±0.03 

3.25 

41.50 ±0.10 

262 

-2.39 

-2.35 ± 0.04 

3.25 

41.70 ±0.10 

128 

-2.71 

-2.61 ± 0.06 

3.25 

41.90 ±0.10 

68 

-2.98 

-2.94 ±0.08 

3.25 

42.10 ±0.10 

28 

-3.37 

-3.17 ±0.13 

3.25 

42.30 ±0.10 

12 

-3.73 

-3.52 ± 0.20 

3.25 

42.50 ±0.10 

3 

-4.34 

-4.12 ± 0.39 

3.25 

z = 1.42 

41.95 ±0.15 

284 

-2.63 

-2.49 ±0.03 

4.06 

42.25 ±0.15 

73 

-3.22 

-3.14 ±0.07 

4.06 

42.55 ±0.15 

12 

-4.01 

-3.89 ±0.19 

4.06 

42.85 ±0.15 

2 

-4.78 

-4.64 ±0.48 

4.06 


z = 2.23 


42.60 ±0.075 
42.75 ±0.075 
42.90 ±0.075 
43.05 ±0.075 

84 

70 

22 

5 

-3.27 

-3.36 

-3.86 

-4.51 

-3.08 ±0.06 
-3.14 ±0.07 
-3.65 ± 0.13 
-4.26 ±0.29 

10.46 

10.69 

10.69 

10.69 

z = 3.24 

42.65 ±0.075 

70 

-3.33 

-3.17 ±0.07 

9.99 

42.80 ±0.075 

52 

-3.48 

-3.26 ±0.09 

10.48 

42.95 ±0.075 

25 

-3.80 

-3.55 ± 0.13 

10.48 

43.10 ±0.075 

6 

-4.42 

-4.17 ±0.27 

10.48 


APPENDIX B: BINNED LUMINOSITY 
FUNCTION 

Here we include two tables that show the binned data points 
of the LF that are plotted in figures and We include 
these plots as a convenience for future studies who wish to 
compare their LFs to ours. 


logio L[oii] 

(erg 

# 

"^obs 

(Mpc-3 d logic L) 

^final 

{Mpc“3 d logic L) 

Volume 
(10^ Mpc^) 

z = 1.47 

41.65 ±0.075 

590 

-2.24 

-2.08 ±0.02 

6.80 

41.80 ±0.075 

425 

-2.38 

-2.28 ±0.03 

6.80 

41.95 ±0.075 

257 

-2.60 

-2.46 ±0.04 

6.80 

42.10 ±0.075 

127 

-2.90 

-2.69 ±0.06 

6.80 

42.25 ± 0.075 

42 

-3.39 

-3.05 ±0.10 

6.80 

42.40 ± 0.075 

19 

-3.73 

-3.55 ±0.15 

6.80 

42.55 ± 0.075 

6 

-4.23 

-4.23 ±0.28 

6.80 

z = 2.25 

42.45 ±0.10 

92 

-3.14 

-2.77 ±0.05 

6.29 

42.65 ±0.10 

37 

-3.53 

-3.15 ±0.08 

6.29 

42.85 ±0.10 

3 

-4.62 

-4.46 ±0.35 

6.29 

z = 3.34 

43.05 ± 0.050 

12 

-4.12 

-3.86 ±0.17 

15.88 

43.15 ±0.075 

7 

-4.37 

-3.92 ±0.24 

16.52 

43.30 ±0.075 

2 

-5.22 

-4.87 ±0.48 

16.52 

z = 4.69 

42.86 ±0.075 

10 

-4.26 

-3.66 ±0.09 

12.22 

43.01 ± 0.075 

5 

-4.56 

-3.93 ±0.13 

12.22 

43.16 ±0.075 

2 

-4.96 

-4.11 ± 0.16 

12.22 


APPENDIX C: STAR-FORMATION RATE 
DENSITY COMPILATION 


In this section, we have compiled a table of the star- 
formation rate densities from different diagnostics spread 
over a wide redshift range. Because each study has its own 
set of assumptions, diagnostics, calibrations, dust correc¬ 
tions,etc. it is quite confusing in keeping track of which study 
has used which set of assumptions. Let alone, for the earli¬ 
est papers, we have to even take into account the different 
cosmologies. To make life much easier for you as the reader 
who may be interested in studying the evolution of the cos¬ 
mic SFR density, we have included in the appendix a long 
table which is our compilation of the SFR densities and lumi¬ 
nosity function parameters from a range of different studies. 
Parts of this table are from Ly et al. (20071, but updated 
with the newest studies in the field. 


This paper has been typeset from a Tj^X/ file prepared 

by the author. 





















24 Khostovan et al. 


Table Cl. SFRD Compilation 


Study 

z 

Diagnostic 


Observed 


logic P* 




loglO 0* 
Mpc “0 

logio 

erg s~^ 

a 

Mq yr“i Mpc“0 


Ciardullo et al. 12013 


Sobral et al.|2012^ 


Bayliss et al^OTl 

Ly et al.|2007 


Zhu, Moustakas fc Blanton|20o'9 


Takahashi et al. 12007 


0.0 - 0.2 
0.2 - 0.325 
0.325 - 0.45 
0.45 - 0.56 
1.47 
1.85 
0.89 
0.91 
1.18 
1.47 
0.84 
1.02 
1.19 
1.37 

1.71 - 1.203 
1.71 - 1.203 


[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 


9 qrv + 0.09 
^•'^'^-0.11 
_ 9 1 9 + 0.05 

^•-L^_0.06 

-2 07+°-°^ 
^•'-'^-0.05 
9 /^rr+0.03 
^•^'-0.08 
- 2.01 ± 0.10 
-2.23 ±0.09 
-2.25 ±0.13 
-1.97 ±0.09 
- 2.20 ± 0.10 
-1.97 ±0.06 


_ 9 qy+O.lO 

^•'^'-0.12 


40 

40.54 ± 0.11 

40.751^0 

40.93in2 

41.71 ± 0.09 
41.31 ± 0.06 
41.33 ±0.09 
41.40 ± 0.07 
41.74 ± 0.07 
41.60 ± 0.05 


41 7Q+0'0'^ 


9 C 7 + 0.28 1. 7C+O.32 


- 1.2 

- 1.2 

- 1.2 

- 1.2 

-0.9 ±0.2 
-1.3 ±0.2 
-1.27± 0.14 
- 1.20 ± 0.10 
-1.15±0.11 
-0.78 ±0.13 


_1 41+0.16 

_ -1 00 + 0.40 

±.00_n 0-7 


Glazebrook et al.|2 

004 

0.90 

[Oil] 

-2.91 

42.30 

-1.3 

'I'epiitz et ai.|2UUo| 

0.90 ±0.50 

[Oil] 

-3.06 ±0.12 

42.15 ±0.08 

-1.35 

Gailego et ai.|2UU2 

0.025 ± 0.025 

[Oil] 

-3.48 ±0.19 

41.24 ± 0.13 

-1.21 ± 0.21 

Hicks et al.|2002| 

1.20 ±0.40 

[Oil] 




Hogg et al.[l998[ 

0.20 ±0.10 

[Oil] 





Hammer et al. 1997 


Colbert et al.|2013 


Pirzkal et al.| 

2013 

Ly et ai.|20U l 



0.40 ±0.10 
0.60 ±0.10 
0.80 ±0.10 
1.00 ±0.10 
1.20 ±0.10 
0.375 ± 0.125 
0.625 ± 0.125 
0.875 ± 0.125 
0.7- 1.5 
1.5 - 2.3 
0.7- 1.5 
1.5 - 2.3 
0.5 ± 0.4 
0.41 
0.42 
0.62 
0.83 


[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[Oil] 

[OIII] 

[OIII] 

[OIII] 

[OIII] 

[OIII] 

[OIII] 

[OIII] 

[OIII] 

[OIII] 


-3.19 ±0.09 
-3.74 ± 0.43 
-3.28 ±0.09 
-3.60 ±0.14 
9 c:o+0.09 
^• 30 - 0.09 

-2.55 ±0.25 
-2.38 ±0.22 
-2.58 ±0.17 
-2.54 ± 0.15 


42.34 ± 0.06 
42.91 ± 0.37 
42.39 ±0.08 
42.83 ±0.11 

41 3+0.09 

^^"^-0.09 
41.17± 0.22 
41.11 ± 0.24 
41.51 ± 0.15 
41.53 ±0.11 


-1.40 ±0.15 
-1.67 ± 0.78 
-1.50 
-1.50 

_ -1 91 +0.08 

^•^^-0.07 
-1.49 ±0.11 
-1.25 ±0.13 
-1.22 ±0.13 
-1.44 ± 0.09 


-2.05 ± 0.11 
-1.82 ± 0.06 
-1.71 ± 0.05 
- 1.66 ±0.06 
-1.48 ±0.10 
-0.92 ± 0.08 
-1.68 ±0.03 
-1.36 ±0.02 
-1.27 ± 0.02 
-1.27 ± 0.02 

— 1 70+011 
1 ycr+O.lS 

l''O-0.08 

— 1 57+0-25 

_ 1 9c:+0.05 

^■^'1-0.08 

— 1 61+000 
1-01-0.28 

_i qc;+0.34 
I'-JO-O.SO 

-1.55 ± 0.06 

-3.02 ± 0.15 

— 1 5Q+0-00 

l-'J7_0.48 

_ 9 07 + 0.11 

^" 1 ' -0.16 

_ 1 77 + 0.09 

1 -' ' - 0.12 

— 1 69”*"®-®® 

l-u»_0.08 

1 7c:+0.07 

l-'O-0.08 

_1 44+0.09 

I'^^-O.ll 

_ 1 C 7 + 0.18 

l'O'-0.30 

—2 20 ’*"®-®^ 
^■^'1-0.08 
1 79+0.11 
l-'^-O.lS 

_ 1 oc:+ 0.20 

i"i'l_o.38 


-2.17 ±0.06 
-2.31 ± 0.09 
-2.06 ±0.05 
-1.73 ±0.03 

























































